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Abstract 
Backdraft Phenomena 
by 
Charles Martin Fleischmann 
Doctor of Philosophy in Fire Protection Engineering Science 


Professor P. J. Pagni, Co-Chair and Professor R. Brady Williamson Co-Chair 


The purpose of this project was to develop a fundamental physical understanding 
of backdraft phenomena. The research was divided into three phases: exploratory 
simulations, gravity current modeling, and quantitative backdraft experiments. The 
primary goal of the first phase was to safely simulate a backdraft in the laboratory. A half- 
residential-scale compartment was built to conduct exploratory experiments. The initial 
experiments concluded with a scenario describing the fundamental physics of backdrafts. 
The importance of the gravity current which enters the compartment after opening was 
identified. In the second phase, the gravity current speed and the extent of its mixed 
region was investigated in a series of scaled salt water experiments. The scaled 
compartment (0.3m x 0.15m x 0.15m) was fitted with a variety of end openings: full, slot, 
door, and window. Video and photo data indicate that the mixing layer which rides on the 
gravity current in the full opening case, expands to occupy nearly the entire current in the 
partial opening cases. The Froude number and nondimensional head height are 
independent of B and are in good agreement with numerical simulations and special limits 
from the literature. 

In the final phase, 28 backdraft experiments were conducted in a 1.2 m by 1.2m 
by 2.4 m compartment. A methane burner was ignited inside a closed compartment and 
allowed to burn as long as oxygen was available. After the flame extinguished due to 
oxygen starvation, the burner was left on to allow the unburned fuel fraction to increase. 


Upon opening the hatch a gravity current enters the compartment and travels across the 
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floor to the ignition source. After ignition a deflagration nps through the compartment 
and out the opening culminating in a large fireball. Histories recorded included: fuel flow 
rates, upper layer temperatures, lower layer temperatures, opening velocities, 
compartment pressures, upper layer species concentrations for O,, CO,, CO, and HC. 
Results indicate that unbumed fuel mass fractions >15% are necessary for a backdraft to 
occur and that the backdraft severity strongly depends on the delay time and species 


concentrations. 


Patrick J. Pagni R. Brady Williamson E: 


Co-Chair, Thesis Committee Co-Chair, Thesis Committee 
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CHAPTER 1 
INTRODUCTION 


1.1 Impetus for this Research 

A backdraft is defined as a rapid deflagration following the introduction of oxygen 
into a compartment filled with accumulated unburned fuel. Prior to 1991, the word 
backdraft was known only to firefighters and a few researchers. However, in 1991 
Universal Studios released a major motion picture called "BACKDRAFT" and almost 
overnight backdraft became a household word. Prior to the movie there was a great deal 
of confusion about the concept of backdraft, and since the movie the situation has become 
even more confused. Articles on backdraft, sometimes called smoke explosions, first 
appeared in the literature in 1914!. This first article described backdraft as a dust 
explosion caused by the carbon particles in the smoke. Other more recently proposed 
"theories explaining backdrafts" require impossible conditions such as instantaneous 
transport of oxidizer to reactant when a compartment is opened, autoignition of gases at 
impossible compartment temperatures, and ignition of soot particles at temperatures less 
than as 500°K.2 

A review of fire service literature reveals many narrative articles on backdrafts and 
the terrible consequences to firefighters caught in the backdraft3:4.°. Typically in the 
articles, firefighters are involved in initial search and rescue or suppression operations 
when the explosion occurs. Although the fire service training manuals specifically provide 
warning signs of backdrafts®.’, in most reported occurrences, the warning signs were not 
observed by the firefighters. Current tactics for reducing the backdraft hazard are to vent 


the structure at the highest possible location prior to entry. However, the ventilation 


process is often a second priority to the rescue operation. Indeed, unless it is restricted to 
roof venting, ventilation may facilitate rather than prevent a backdraft. Lack of warning 
and poor ventilation practices are commonly reported in descriptions of backdraft. 

As a result of a backdraft which killed two firefighters in Chatham, England in 
19753, an extensive literature review on fires involving explosions was published by 
Croft®. This review covered the United Kingdom, the United States, and Canada during a 
seventy year period ending in 1976. A total of 127 fires involving explosions were 
reported in the literature with sufficient detail to be included in this study. The explosions 
were categorized as occurring in smoldering fires (52), developing fires (57), and 
developing fires with secondary explosions (18). The latter category is not relevant to 


backdrafts. Important conclusions from Croft's work are?: 


1) More explosions were associated with developing fires than with 
smoldering fires. 


2) More fire fighting personnel were injured in explosions associated with 
developing fires, than with smoldering fires. 


3) A large number of fire fighting personnel were killed in explosions 
associated with smoldering fires because the firefighters had entered the 
premises believing that the fire had been quelled or that it had extinguished 
itself due to a lack of oxygen. 


4) More fire fighting personnel were injured in explosions in shops and 
supermarkets than in any other type of occupancy. 


5) More fire fighting personnel were killed in explosions in warehouses 
than in any other single occupancy. 


6) Smoke explosions are more likely to occur in factories and warehouses 
than in any other occupancies. 


7) More explosions resulted from fires involving cellulosic materials than 
any other material. This is most probably associated with the pattern of 
usage of different materials during the period studied. 


8) The cellulosic materials which caused the largest number of explosions were 
varnished, painted, and polished woodwork and combustible fiberboard. 


The results presented in Croft's study may be distorted. Although his study was extensive 
it was not comprehensive because it was limited to the available literature. Accurate 
statistical data were not available. 

National statistics on the number of firefighters killed or injured as a result of 
backdrafts are not now available. The fire analysis and research division of the National 
Fire Protection Association does not identify backdraft as a separate cause of firefighters’ 
death or injury; it simply combines backdraft casualties with others under the heading of 
"Rapid Fire Progress in Structures". The reason for this lack of distinction is that it is 
often difficult to identify the specific phenomena actually occurring during a fire’. Other 
phenomena in this category are flashover and ignition of fire gases. Other statistics on 
backdrafts such as the annual number of occurrences and number of civilian casualties are 
also not available. 

In addition to the literature review, the Chatham fire also prompted some 
experiments. Smoldering foamed rubber was placed in a small enclosure, 1.4 m3, and 
allowed to accumulate excess pyrolyzates!°. When a diffusion flame was introduced near 
the bottom of the enclosure, an explosion occurred!!. This was the only research on 


backdraft that this author could find within the commonly available literature. 


1.2 Backdraft Scenario 

As a result of the exploratory experiments presented in Chapter 2, a fundamentally 
sound backdraft scenario has been identified. Consider a fire in a closed compartment 
where only a minimal amount ventilation is provided by leakage. As the fire heats the 
compartment, leaks in the compartment bounding surfaces permit outflows that minimize 
any pressure differential!2, A hot layer composed largely of combustion products 


descends around the fire causing some pyrolysis products to remain unbumed. These 


products accumulate forming a deep layer which has insufficient oxygen to support 
combustion and becomes fuel rich as the fire continues to supply unburned fuel. We 
assume a small flame or glowing ember remains burning somewhere in the lower layer. 
Suddenly, a new ventilation source is provided by a window breaking}? or door opening. 
The hot, vitiated atmosphere within the compartment flows out the upper portion of the 
vent. Simultaneously, cold, fresh air flows in the lower portion of the vent. This cold, 
density driven, flow is called a gravity current!*. A mixed layer forms due to the 
instabilities at the shear interface between the outflow and the inflow and rides on the 
gravity current across the compartment!5. For many opening geometeries, large scale 
vortical structures with fuel rich and oxygen rich regions rolled up sequentially, fill the 
gravity current. Figure 1.1a shows a salt water model of a gravity current nearly 
completing its propagation across a compartment. A portion of this mixed current is 
within the flammable range and is ignited when it reaches a flame or glowing ember. After 
ignition, a new flame propagates back through the mixed region. Figure 1.1 b shows a 
flame approximately 1 s after ignition in an experimental compartment. The flame shape 
shows a precursor flame leading the primary burning zone. This preceding flame burns 
along the interface between the entering gravity current and the exiting compartment 
gases. The preceding flame and its wake are sufficiently unstable to generate a rapidly 
propagating turbulent flame. The resulting turbulent deflagration within the compartment 
drives some accumulated excess pyrolyzates out the vent and consumes that fuel outside 
the compartment in the dramatic fireball commonly associated with backdraft. This entire 
process: the accumulation of unburned gaseous fuel, the propagation of an oxygen rich 
gravity current creating a mixed region and carrying it to the ignition source, the ignition 
and propagation of an eventually turbulent deflagration and the external fireball, altogether 


constitutes a backdraft. 
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Figure 1.1a - Photograph of the salt water model showing the entering gravity current. 
B=0.080, opening was the h,/3 horizontal slot, ~2.4 s after opening. 


Figure 1.1b - Photograph of the flame propagation along the top of the entering gravity 
current. The ignition spark is turned when the compartment is opened. Compartment is 
1.2 m wide by 2.4 m long by 1.2 m high and the opening is 1.1 m wide by 0.4 m high. 


Two modifications to the above scenario were observed. Identical compartment 
conditions at opening are assumed but there is an increase in the time delay before ignition 
of the backdraft. This time delay results in a different deflagration flame structure and 
backdraft severity. In the first modification, the time delay is sufficiently long that the 
gravity current may have even returned to the front opening and the compartment will act 
like a reservoir in a reservoir filling problem. Figure 1.2a shows the current as it returns to 
the opening in the salt water model. The lower layer formed by the gravity current is 
made up of entrained fuel rich compartment gases along with the oxygen rich air entering 
the compartment. If ignition occurs during this time period, with the lower layer more 
uniformly mixed and within a flammable range, the flame structure is hemispherical in 
shape and the initial deflagration is more severe. Figure 1.2b shows the hemispherical 
flame identifying this scenario. 

The second modification requires a long ignition delay where the gravity current 
has left the compartment and the lower layer is made up primarily of air. The flammable 
compartment gases are trapped above the soffit and the lower layer is primarily made up 
of air. Figure 1.3a shows the salt water model of the compartment with the gases trapped 
above the soffit. If ignition occurs now, a flame will travel along the interface as seen in a 
mine gallery!®. The propagating flame and its wake are sufficiently unstable to generate a 
rapidly propagating turbulent, albeit less severe, deflagration. Figure 1.3b shows the flame 
propagating along the interface. In this experiment ignition occurred 300 s after vent 


opening. 
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Figure 1.2a - Photograph of the salt water model showing the reflected gravity current. 
B=0.080, opening was the h,/3 horizontal slot, ~6.0 s after opening. 
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Figure 1.2b - Photograph of the hemispherical flame which is created when ignition is 
delayed until the gravity current was reflected off the wall opposite the opening. 


Compartment is 1.2 m wide by 2.4 m long by 1.2 m high and opening 1s 1.1 m wide by 
0.4 m high. 
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Figure 1.3a - Photograph of salt water model showing the compartment fluid trapped 
above the soffit. B=0.080, opening was the h,/3 horizontal slot, ~60 s after opening. 


Figure 1.3b - Photograph showing the flame propagating along the interface between the 
lower layer, made up primarily of air, and the fuel rich upper layer trapped above the 
soffit. Flame was ignited by a spark 300 s after opening. The spark ignitor was located 
0.15 m from the wall opposite the opening and at the height of the soffit. Compartment is 
1.2 m wide by 2.4 m long by 1.2 m high and the opening is 1.1 m wide by 0.4 m high. 


1.3 Research Outline 

Since little research on backdraft is available, the first step was to formulate a 
working scenario as a physical explanation of the phenomena. A series of exploratory 
experiments were conducted in a special compartment designed to safely simulate 
backdraft experiments. A single horizontal slot opening was investigated. Description of 
the apparatus and results of these experiments are presented in Chapter 2 along with 
computer modeling results which were used to investigate the species concentrations 
within the compartment prior to backdraft. 

Once the feasibility of the scenario was demonstrated, a series of salt water 
experiments were conducted to study the geometry, mixing, and velocity of the gravity 
current which enters the compartment prior to backdraft. The results of the salt water 
experiments are given in Chapter 3. The velocity and height of the entering gravity 
current are presented. In addition, a two dimensional numerical simulations of the gravity 
current experiments were performed at the National Institute for Standards and 
Technology. These numerical predictions are compared with the salt water modeling and 
backdraft experimental results in Chapter 4. 

Chapter 5 discusses results from a second series of backdrafts which were more 
fully instrumented than the exploratory experiments. Data collected include burner flow 
rates, compartment temperatures, species concentrations, layer height, vent flow rates, and 
compartment pressure. Exemplar data are presented along with a summary table for the 
experiments using a horizontal slot opening. A series of experiments were also conducted 
using a simulated window opening. The results of these experiments are given in Chapter 
6. The fractions of the chemical energy stored as unburned hydrocarbons within the 
compartment at the opening time which were released inside and outside the compartment 
are also discussed in this chapter. Chapter 7 gives the conclusions of this research and 


discusses possible directions for future work. Appendix A gives the details of the species 


balance performed on the gas samples to determine the concentrations of water and 
nitrogen and as a check on the hydrocarbon concentration. Details of the experimental 
apparatus are covered in appendix B including design drawings and construction material 
information sheets. Appendix C is a summary of the data collected in all of the backdraft 
experiments. Although all of the experiments performed are not discussed in this work, all 


of experiments are included here for completeness and future reference. 
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CHAPTER 2 
EXPLORATORY BACKDRAFT EXPERIMENTS 


2.1 Introduction 

Fires can produce more fuel than the locally available oxygen can consume. This 
surplus fuel is called excess pyrolyzates!. If the compartment containing the fire is well- 
ventilated, the excess pyrolyzates fuel long flames that extend out openings in the 
compartment, rapidly spreading the fire?. If the compartment is closed, the excess 
pyrolyzates accumulate, ready to burn when a vent is suddenly opened, e.g., by a window 
breaking due to the fire-induced thermal stress? or by a firefighter entering the 
compartment>-*78, Upon venting, a gravity current carries fresh air into the compartment. 
This air mixes with the excess pyrolyzates producing a flammable pre-mixed gas which 
can be ignited in many ways. The rapid deflagration moving through the compartment 
after ignition, consuming the accumulated excess pyrolyzates, is called a backdraft. 

The fire service community has long recognized the hazards associated with 
backdrafts?!°1!_ The literature provides a definition of backdraft!2: "Backdraft is the 
burning of heated gaseous products of combustion when oxygen is introduced into an 
environment that has a depleted supply of oxygen due to fire. This burning often occurs 
with explosive force." This definition would be correct if the word pyrolysis were 
substituted for combustion since it is primarily the unburned pyrolyzates which provide the 
fuel for backdrafts. 

None of the quantitative models of compartment fires currently available!3-!4 15 
incorporate backdraft phenomena because the underlying fundamentals are not well 


understood. This research is aimed at advancing our understanding of the physics and 
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chemistry of backdrafts. A backdraft scenario is presented. Half-scale room fire 
experiments attempting to confirm this scenario and results of computer modeling of pre- 


backdraft compartment conditions are described. 


2.2 Backdraft Scenario 

Consider a fire in a closed compartment where the only ventilation provided is by 
leakage. As the fire heats the compartment, leaks in the compartment bounding surfaces 
permit outflows that minimize any pressure differential!®. A hot layer composed largely of 
combustion products descends around the fire causing some fuel pyrolysis products to 
remain unburned. These products accumulate forming a deep, fuel-rich layer. We assume a 
small flame or glowing ember remains burning. Suddenly a new vent is opened. The hot, 
vitiated atmosphere within the compartment flows out of the upper portion of the vent. 
Simultaneously, cold, fresh air flows into the lower portion of the vent. The propagation 
of the leading edge of this cold, density driven, flow is called a gravity current!’. Figure 2- 
1 shows a gravity current flowing into a compartment in a salt water model. A mixed 
layer forms due to instabilities at the shear interfaces between the outflow and the inflow, 
and is carried across the compartment by the gravity current!®!9_ A portion of this mixed 
current is within the flammable range and is ignited when it reaches a flame or glowing 
ember. After ignition, the new flame propagates back through the mixed layer. That flame 
and its wake are sufficiently unstable to generate a rapidly propagating turbulent flame. 
The resulting turbulent deflagration within the compartment drives any remaining 
unburned fuel and pyrolysis products out the opening to burn outside the compartment in 


a dramatic fireball. 


Figure a We Photograph of a gravity current in the salt water modeling experiments 
replicating the backdraft compartment as described in chapter 3. B = 0.023, h,/3 
horizontal slot opening, ~5 s after opening compartment. 


2.3 Description of Experimental Apparatus 

To test the hypothesized physical explanation of backdraft an experimental 
program was undertaken. The primary goal of this program was to safely simulate 
backdrafts in the laboratory. 

Because of the explosive nature of backdraft, a special chamber was constructed to 
replicate a small room at approximately half scale. Figure 2-2 shows a schematic of the 
apparatus giving the internal dimensions of the compartment. Figure 2-3 is a photo of the 
apparatus. In order to control the overpressure hazard, one long wall was designed as a 
pressure relief panel. The entire wall was hinged along the bottom and closed with a single 
nylon fastener at the top. Failure of the fastener was designed to relieve any overpressure 
greater than 1 kPa, as recommended for venting?®2!. Tests with a large, pressurized 


plastic bag showed that the blow-out panel released at 0.6 + 0.1 kPa. It is also 


14 


recommended that the pressure relief panel weight be limited to <15 kg/m? to reduce 
inertia and opening time. The pressure relief wall weighed approximately 13 kg/m. It was 
constructed of 1.2 mm (18 gauge) steel studs, 50 mm wide, 1.6 m long, 0.6 m on center. 
The sheathing over the studs was 1.2 mm galvanized sheet steel. The panel interior was 
covered with a 25 mm thick layer of refractory fiber blanket. 

All the stationary walls, ceiling, and floor were designed to withstand 5 kPa, five 
times the expected maximum overpressure. The three stationary walls used 1.2 mm steel 
studs, 0.10 m wide, 1.4 m long, 0.4 m on center. The ceiling was constructed of 1.2 mm 
steel studs, 0.2 m wide, 1.4 m long, 0.4 m on center. Two layers of 16 mm Type X, fire 
rated, gypsum wallboard were mounted as interior sheathing to the stationary walls and 
ceiling to provide structural strength and secondary thermal fire resistance. Gypsum 
wallboard was also installed on the floor to protect the plywood platform. A 50 mm thick 
refractory fiber blanket was installed over the gypsum wallboard on ihe walls and ceiling 
to provide the primary thermal resistance for the structure. This insulation allows for 
repeated experiments without the need to rebuild the chamber. 

A 0.9 m high and 1.5 m wide observation window was installed in the wall 
opposite the pressure relief panels, as shown in Fig. 2-3. The window glass was 
Neoceram??, a ceramic with a negative coefficient of expansion below 900K and is capable 
of withstanding continuous exposure to temperatures of 1000K. The glass was mounted in 
a standard steel frame protected from the hot compartment gases by a refractory insulation 


blanket. 
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Pressure Relief 
Panel 


Opening Hatch qi 


Figure 2.2 - A schematic diagram of the half- room-scale backdraft apparatus showing 
dimensions and component locations. 


Figure 2.3 - Photograph showing the observation window in the backdraft apparatus. 
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To simulate a window or door, a 0.4 m high, 1.2 m wide opening was centered in 
the short wall opposite the burner, see Figs. 2-2 and 2-3. This vent was covered with a 
Both the pressure relief panel and the opening hatch are in the open position. 
manually operated hatch which was opened after the fire had been burning for several 
minutes. The hatch was hinged at the bottom and held closed by a single throw latch at the 
top. Figures 2-2 and 2-3 show the hatch in the open position. 

A gas burner, 0.3 m square and 0.3 m high, was used in all these experiments. A 
spark ignitor mounted 50 mm above the burner, centered on the edge toward the 
compartment center, was the ignition source for both the burner and the backdraft. A 
10,000 volt transformer was used to produce the arc between two electrodes 5 mm apart. 
The burner was placed against the wall opposite the opening, as seen in Fig. 2:9) Every 
effort was made to seal all construction holes to control leakage. Two small rectangular 
vents, 25 mm high, 0.3 m wide, were placed at the floor and ceiling in the wall with the 
hatch to allow for controlled leakage. In most experiments only the floor vent was open. A 
vertical thermocouple tree was placed in the geometric center of the compartment, as 
shown in Fig. 2-2. The thermocouples were 0.5 mm type K thermocouple wire with a 
stainless steel overbraid. The average bead diameter was 1.1 mm. Seven thermocouples 
were located at 0.15 m intervals, measured from the ceiling. An additional thermocouple 


was placed 50 mm from the ceiling to measure the ceiling jet temperature. 


2.4 Results 

A total of 23 backdraft experiments were conducted. The experimental parameters 
are summarized in columns 2 through 6 of Table 2-1 including: fuel, fuel flow rate, burner 
time, opening time, and number of leakage vents. The experiment numbers given in 
column | are used throughout this paper. The last column in Table 2-1 indicates if a large 


fire ball was observed outside the compartment. Density relative to air was the criterion 
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for selecting the two gaseous fuels, natural gas and propane. Different fuel flow rates, 
burn times, and ventilation configurations were used in order to vary the excess 
pyrolyzates (unburned fuel) stored within the upper layer. Two different ventilation 
configurations were used. In experiments 1 through 11 and 20 through 23 a single vent at 
the floor provided leakage. In experiments 12 through 19, an additional vent was placed 
at the ceiling. None of the two vent experiments resulted in a backdraft. 

Only 8 of the 23 experiments resulted in backdrafts, i.e., 1 - 8. The experiments 
were not considered to result in a backdraft unless a large fire ball was observed outside 
the compartment. The 8 experiments which resulted in a backdraft used natural gas as 
fuel. In experiments 4, 5 and 7, the experimental parameters were constant as indicated in 
Table 2-1: a burner time of 175 s, a single floor vent and a 5 s delay between burner shut 
off and the opening of the hatch to allow the fluid mechanics caused by the burner to 
subside. In experiment 1 the burn time was shorter at 170 s and a2 s delay. The burn 
times in experiments 6 and 8 were 180 s and 185 s, respectively. 

In experiment 9, a similar flame structure was seen but it did not result in a large 
fire ball. When opening the compartment, the latch was released and the hatch opened 
slightly but did not fall open. The hatch was then pulled open by the operator. This 
manner of opening allowed flow into the compartment ahead of the gravity current. 
Ignition had occurred before the incoming gravity current reached the ignitor, causing the 
compartment pressure to increase and thus reduce the incoming air flow and lowering the 
energy release. 

In experiments 10 and 11 the burner was on at the time of opening which caused 
additional mixing in the area of the ignitor and interfered with the gravity current ignition. 
When the ceiling vent was opened as in experiment 12, insufficient hydrocarbons were 


available to fuel the backdraft. 


18 


Table 2.1 - Experimental parameters and results from exploratory backdraft experiments. 


Fuel Burner | Opening 
Experiment Flow rate Time Time Number Fire 
Samet Fuel (kW) (s) (s) of vents Ball 
yes 


Natural Gas 
Natural Gas 
Natural Gas 
Natural Gas 
Natural Gas 
Natural Gas 
Natural Gas 
Natural Gas 
Natural Gas 
Natural Gas 
Natural Gas 
Natural Gas 
Propane 
Propane 
Propane 
Propane 
Propane 
Propane 
Propane 
Propane 
Propane 
Propane 
Propane 
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* Glare on the observation window washed video recording 


In the propane experiments, experiments 13 through 23, the density of the propane 
relative to the other compartment gases resulted in the accumulation of the propane low in 
the compartment. With the fuel near the floor, a small dancing flame stabilized near the 


small open floor vent. When the hatch was opened, the fuel was quickly ignited by the 
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dancing flame and burned as dark orange and yellow flames which slowly propagated 
throughout. 

The data collected in these experiments has been limited to the temperatures 
measured on the thermocouple tree and data recorded through the window using 35 mm 
cameras and video camcorders. Typical temperature histories at different heights within 
the compartment for experiment 4 are shown in Fig. 2-4. The burner was ignited at time 
zero. The temperature at the top of the compartment rose quickly to a maximum of 820 K 
at 25 seconds after ignition. After this peak the temperature dropped as the compartment 
lost energy through its boundary surfaces and the burning rate was limited by the lack of 
oxygen. After 120 seconds when the fire was nearly out the flames detached from the 
burner and began to dance across the floor, as seen in Fig. 2-5a. The dancing lasted ~ 30 
seconds and was responsible for the temperature rise shown between 120 and 150 seconds 
in Fig. 2-4. These dancing flames occurred in some but not all of the experiments. Similar 
behavior has been described by Sugawa et al.?3 in their work on poorly ventilated pool 
fires within a compartment. Figure 2-5b, taken ~ 3.5 s after the vent is opened, just after 
the gravity current reaches the spark, shows the propagation of a mostly premixed flame 
through the mixed layer formed between the hot, fuel-rich, upper layer gases and the cold, 
oxygen-rich, fresh air entering the compartment through the lower portion of the open 
vent. Similar premixed flames have been reported by Phillips?* on a buoyant natural gas 
layer interface within a model mine gallery, open at the bottom to allow free expansion of 
the combustion products. Phillips identified three flames: a premixed U-shaped flame 
burning where flammable natural gas concentrations occurred, a diffusion flame at the 
natural gas/air interface behind the premixed flame, and an unstable flame formed in the 
hot product layer sandwiched between the cold natural gas and air layers. In the backdraft 
experiments described in this paper, the burning occurs within a closed chamber which 


restrains the hot products. As the burnt gases expand, they force the unburned fuel and air 
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ahead of the advancing flame front, out the vent. This behavior is demonstrated by the 
large fire ball which burns outside the compartment, shown in Fig. 2-5c. The spike in the 
temperature, seen in Fig. 2-4 after 180 s, is the deflagration wave flame front as it moves 


past the thermocouple tree on its way out of the compartment. 


950 
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Figure 2.4 - Representative compartment temperature histories for Experiment 4 in Table 
2.1. Locations are measured from the floor ( -) 1.02 m, (-——) 0.72 m, (------) 0.42 m, 
and (——-) 0.12 m. 
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Figure 2.5a - Photograph showing the dancing flame ~130 s after ignition of the burner in 
Experiment 4 in Table 2.1. 


Figure 2.5b - Photograph showing the premixed flame in the mixed region at the interface 
between hot fuel rich and cold oxygen rich layers in Experiment 5 in Table 2.1. 
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Figure 2.Sc - Photograph showing the large fireball bursting out of the compartment in 
Experiment 8 in Table 2.1. 


In Table 2-2, the ambient temperature, average compartment temperature, 
calculated gravity current travel time, ignition delay time, deflagration wave travel time, 
and total travel time are shown in columns 2 through 7, respectively. The experiment 
numbers given in column 1 correspond to those in Table 2-1. In experiments 2 and 3 the 
sun's glare on the window washed out the video camera image. The ignition delay time 
and deflagration travel time are determined from the video tapes by counting the individual 
frames at 30 frames per second. Taking into account all errors, accuracy is conservatively 
estimated at + 0.2 s. The time of opening was taken as the time when the hatch was at 45° 
with the horizontal. This was done to avoid including slight variations in the hatch 
movement upon opening. The time from the opening of the compartment to ignition, 
shown in column 5, ranges from 2.1 s to 5.2 s. The gravity current travel time, shown in 


column 4, is calculated from gravity current velocity, 
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] 
Ne Beh, (2-1) 


where v is the speed of the leading edge of the gravity current, B = Ap/p, with Ap is the 
initial density difference across the opening, p is the compartment density, h is the 
compartment height, and g is 9.8 m/s?. This expression is developed in chapter 3. 
Densities were calculated using the average compartment temperature at the opening time 
and the ambient temperature. 

The deflagration wave travel time, shown in column 6, is the time from ignition to 
the time the leading edge of the flame leaves the compartment; it ranges from 1.2 s to 2.6 
s. In experiment 8, the ignition delay time was significantly longer, and the premixed flame 
was not observed. When ignition finally occurred, the flames immediately filled a large 
turbulent hemisphere. It should also be noted that the deflagration travel time was 
considerably faster than for previous experiments. The reason for the delayed ignition 
may be that the initial gravity current was still too rich and additional flow into the 
compartment was required to produce a mixed layer within the flammable range. The last 
column in Table 2-2 is the total time required for the gravity current to move across the 


compartment, ignition to occur, and the deflagration to leave the compartment. 


Table 2.2 - Travel time determined from videotapes and gravity current velocities. 


Average Gravity | Ignition} Deflagration | Total 
Ambient Compartment 
Experiment | Temperature | Temperature 
Number 
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The variation of the ignition times, shown in Table 2-2, for similar experimental 
conditions may be caused by nonuniformities within the gravity current head. Calculations 
discussed in chapter 4 show highly non-uniform large scale vortical structures within the 
gravity current. While the mixed region may be flammable on average, it may not ignite 


because the fuel and oxygen have not yet mixed. 


2.5 Backdraft Compartment Fire Modeling 

The computer program FIRe Simulation Technique, FIRST?°, was used to analyze 
the compartment conditions by predicting the histories of the upper layer temperature, 
depth, and species concentrations. The FIRST model is a direct descendant of the Harvard 
compartment fire code?6 and is currently available from the National Institute of Standards 
and Technology in Gaithersburg, Maryland, USA. It provides a time dependent solution of 
simplified species and energy equations governing a compartment fire. The model assumes 
that the compartment can be broken into four large zones, the upper layer, the lower layer, 
the fire plume and the burning objects. Detailed explanations of compartment fire zone 
modeling are available in the literature?’. 

The input data were well defined since all difficult parameters such as leakage rate 
and rate of heat release were specified in the experiments. The following assumptions 
were made: 1) The fire source was a gas burner, 0.34 m in diameter, flowing natural gas at 
150 kW, with an initial 100% efficiency, placed in a corner. FIRST requires a circular 
burner, therefore an equivalent diameter was used. The corner location was chosen since 
the plume entrainment model used in FIRST overestimates entrainment close to burner”. 
2) The only leakage into and out of the compartment was through the 2.5 cm high, 30 cm 
wide vent at the floor. 3) The thermal response of the compartment was governed by the 
refractory fiber blanket. 4) A list of all the other assumptions inherent to FIRST is 
available??. 


25 


Figure 2-6 shows a comparison of the FIRST upper layer temperature history with 
data from experiment 5. The experimental temperature histories from the thermocouple 
tree, e.g. see Fig. 2-4, were converted into unsteady average upper and lower layer 
temperatures using the method Quintiere et. al.3° applied to steady state temperature 
profiles. Quintiere assumed an upper layer temperature from the data and then used the 


following equations to solve for the lower layer temperature and the thermal layer 


interface height: 
h, 1 
{(= }oe= —hy]/T™ +h, /T, (2.2) 
0 
h, 
and = | Tdx =[h,-h, JT’ +h, T”, (2.3) 


0 


where TUL and T/! are the upper and lower layer temperatures, and h, and h, are the 
heights of the compartment and the layer interface. Equation (2-2) is a mass balance and 
Eq. (2-3) retains the same mean temperature as in the data. Here, the lower layer 
temperature was specified as the arithmetic average of the two lower thermocouples and 
the upper layer temperature and thermal interface location were calculated from Eqs. (2-2 
and 2-3). Comparisons for experiment 5 are shown in Figs. 2-6a and b at 60 s and 180s 
respectively. The upper layer temperature and thermal interface height data are indicated 
by the x's in Fig. 2-7a & b, respectively. Both the temperature and thermal interface 
compare well with the FIRST results shown as solid lines. The oxygen and unburned 
hydrocarbon mass fraction histories calculated by FIRST are shown in Fig. 2-7c. Future 
experiments will obtain species concentrations in the upper layer. The hatched vent is 
opened at 180 s in the modeling, causing the rise in oxygen and decrease in hydrocarbons 


shown in Fig 2-7c. The vitiated layer also rises and cools, as shown in Figs. 2-7b and 2-7a, 
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calculations are only applicable to the point of ignition of the backdraft. The experimental 


temperature spikes shown in Figs. 2-4 and 2-7a are not depicted by the computer model. 


Height (m) 


300 350 400 450 500 550 600 650 
Temperature (K) 


Figure 2.6a - Comparison of temperature data from Experiment 5 (x) (Table 2.1) at 60s 
with the idealized two layer approximation calculated from Equation (2-2) and (2-3). 
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Figure 2.6b Comparison of temperature data from Experiment 5 (x) (Table 2.1) at 180 s 
with the idealized two layer approximation calculated from Equations (2-2) and (2-3). 
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Figure 2.7a - Comparison of temperature data (x) from Experiment 5 (Table 2.1) 
calculated from Equations (2-2) and (2-3) with computer results from FIRST (——). 
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Figure 2-7b - Comparison of Experiment 5 (Table 2.1) data (x), calculated from Equations 
(2-2) with the thermal interface history computed from FIRST ). 


28 


Mass Fraction (kg/kg) 


Time (s) 


Figure 2.7c - Compartment upper layer hydrocarbon (— — —) and oxygen ( 
fraction computed by FIRST. 


mass 


2.6 Discussion 

Each aspect of the hypothesized backdraft scenario appears to have been observed 
in these experiments. Before the compartment was opened there were high concentrations 
of unburned hydrocarbons in the upper layer, the temperature was relatively cool and the 
upper layer nearly filled the compartment. After the hatch was opened there was a delay 
of at least 2 seconds before ignition. This delay is the time required for the gravity current 
to travel to the ignition source. Once ignited, a small premixed flame front travels into the 
flammable mixed layer preceding the large non-premixed deflagration. Figure 2-8 is a 
series of video images showing the backdraft development from ignition for experiment 5 
in Tables 2-1 & 2-2. The images are taken at 0.4 s intervals. In the first three images, the 
flame in the mixed region can be clearly seen. Behind this flame is a turbulent combustion 


region which develops from the buoyant instability of the hot combustion products, as 
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shown in the next three images. The hot products rise and displace the fuel-rich upper 
layer downward into the oxygen-rich air. Combustion then further enhances the mixing. 
As the quasi-premixed flame approaches the center of the compartment, Fig. 2-8, the 
expanding turbulent flame has accelerated sufficiently to engulf the laminar flame and 
produce a deflagration which advances rapidly through the compartment and bursts out of 


the opening in a dramatic ball of flame, see in image six. 
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2.7 Conclusions 

The physical model postulated here appears to accurately describe the backdraft 
phenomena observed in these half-scale experiments. Ignition does not occur immediately 
upon opening a vent to a fuel rich compartment. For a fuel rich backdraft to occur, a 
gravity current, on which a mixed layer rides to an ignition source, is required. It is this 
time delay caused by the gravity current propagation which creates a hazard to firefighters 
who may enter a compartment and become trapped in the backdraft process. 

Future work will focus on improved instrumentation in the compartment. A 
variety of opening geometries will be examined. Salt water modeling has been performed 
to determine the size and location of the mixing layer as a function of the density 
difference and the opening configuration. Future compartment experiments will also be 
more fully instrumented with additional thermocouples, pressure transducers, analysis of 
the upper layer gas concentrations including HC, O07, CO, and CO, bi-directional velocity 
probes in the hatch and vent openings, and improved video recording for better flow 
visualization. Hopefully, a foundation will be laid for more sophisticated compartment fire 


models which can incorporate backdraft phenomena. 
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CHAPTER 3 


SALT WATER MODELING OF FIRE COMPARTMENT GRAVITY CURRENTS 


3.1 Introduction 

A gravity current is the flow of one fluid into another caused by a difference in 
density. This density difference may be due to a dissolved chemical or a difference in the 
temperature between the two fluids. There are many common examples of gravity 
currents including sea-breeze fronts, avalanches, lock exchanges, and flows following 
volcanic eruptions. A large body of research is available on the subject of gravity 
currents.!-2_ In many cases, the flow field in gravity currents is sufficiently complex that 
the problem is difficult to solve from first principles. For this reason, physical models, 
typically salt water models, are used to analyze these problems. Salt water models have 
been applied to many fire problems including corridor smoke flow, ship board fires, and 
compartment fires.3>4° 

When a fire occurs in a closed compartment where the only ventilation is due to 
leakage, the fire initially burns independent of the surroundings and a hot upper layer 
develops within the compartment. If the leakage rate is small, the hot layer descends over 
the fire and the burning becomes limited by the available oxygen thus producing large 
amounts of unburned fuel. Left undisturbed, the heat release rate will decrease and the 
fire may enter a smoldering stage. When the compartment is opened, a gravity current of 
dense ambient air enters the compartment mixing with the lighter, fuel rich, compartment 
gases. If the fuel concentrations are high enough the mixed region carried with the gravity 


current may ignite, resulting in a backdraft. The important role of gravity currents in 
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backdrafts has been demonstrated in a series of half-scale experiments discussed in chapter 
Hy 

This study attempts to quantify the gravity current which enters a compartment 
prior to a backdraft. A salt water scale model with two different density fluids is used to 
visualize the flow into the compartment. Because the fire is assumed to be small or 
smoldering, plume effects are ignored and the compartment is filled with a uniform density 
fluid lighter than the fluid outside the compartment. Data collected from these 
experiments include gravity current propagation velocity and geometry. In addition to the 
entering current, the current which is.reflected off the wall opposite to the opening wall is 


examined. 


3.2 Gravity Current Scaling 

As a simplified limit of a compartment fire gravity current, consider the steady 
flow of a perfect fluid in a semi-infinite horizontal box of arbitrary width as shown in Fig. 
3.1. At time zero, far to the nght, the entire end of the box is instantaneously removed. 
High density ambient fluid, state 0, flows into the box, as low density compartment fluid, 
states 1 and 2, flows out, due to buoyancy. The parameter indexing that buoyancy is the 


normalized positive density difference, 


_ (Po =P) 
B aie (3.1) 


where py, is the higher density, ambient, fluid within the gravity current and p , 1s the lower 


density fluid, (at opening) within the hot compartment ahead of the gravity current. 
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Figure 3.1 - Gravity current schematic. Velocities are indicated in a reference frame fixed 
on the gravity current. Heights are indicated by h. 


Benjamin! has shown that, in this perfect fluid limit with no mixing or dissipation, 


conservation of mass, momentum and energy can be written , respectively, as 


v,h, =v,h,, (3.2) 
v, h, +Bgh,° = 2v,"h, +Bgh,”, (3.3) 
v, = 2Bg(h, —h,). (3.4) 


Eliminating v, from Eqs. 2 and 3 and v, from Eqs. 3 and 4, the height of the exiting 


compartment fluid is 


h 
hy = = hy. (3.5) 


This is also the height of the gravity current since h) = h,-h,. Benjamin! shows that for 
flows with energy losses, h)<h,/2. The nondimensional velocity or Froude number of the 


fluid exiting the compartment, from Eqs. 4 and 5 is 


V, — 
ire BS 


since this is >1 a dissipative hydraulic jump is possible. The velocity of the compartment 
fluid approaching the gravity current, or in the laboratory reference frame, the gravity 


current velocity, from Eqs. 2, 5 and 6, is, 
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peeves (3.7) 
Bgh, 2 
Thus for scaling, the characteristic dimension and velocity are 
x = heandive—-)/pehe (3.8) 
from Eqs. 5 and 7 respectively. The characteristic time is then 
(Xen ity (3.9) 
Ve eS 


Typical dwellings have room heights of 2.4 m (8 ft), so the compartment used for 
backdraft experiments, height 1.2 m (4 ft), is 1/2 scale. The salt water compartment 
described here is 0.15 m (0.5 ft) high, so it is 1/8 scale to the experimental backdraft 
compartment and 1/16 scale to a dwelling. 

The salt water experiments are necessary to quantify the effect of transients, 
mixing, energy dissipation, opening geometry, and aspect ratio on the simple gravity 
current size and nondimensional velocity expressions given by Eqs. 5 and 7. They will 
provide confirmation and corrected formulas useful for modeling backdrafts in fire 
compartments. | 

The salt water experiments are limited to 0.003 < B < 0.101, while the backdraft 
compartment and full scale fires produce higher B, up to 1.2. However, the literature 
suggests,® as confirmed by experimental Froude numbers developed here, that v* is 
independent of B and depends only on the opening geometry. Therefore, these eealed 
velocity and geometry results are expected to apply directly to actual and modeled 
backdrafts. For example, a velocity of 0.09 m/s at a B = 0.05 for a salt water current in a 
slot opening geometry gives v* = 0.32, which would correspond to 0.8 m/s at a B = 0.5 in 


the 1.2 m high model compartment and 1.1 m/.s at the same B in a 2.4 m high dwelling. 
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Heat transfer effects at the boundaries and between fluids are not included‘. 
Separate analyses and experiments will be required for gravity currents submerged at great 
depths |? as would occur upon opening a small door in a large, high warehouse with a 


ventilation limited fire. 


3.3 Apparatus and Procedure 

Salt water experiments were conducted by placing an acrylic compartment within a 
larger glass tank. The large tank, 0.3 m wide, 0.6 m long, and 0.45 m deep, contained a 
dense water and salt solution ranging in density from 1.003 kg/m? to 1.101 kg/m3. The 
solution temperature was 18°C. Standard "Rock Salt" crystals were dissolved in tap water 
to raise the density to the desired level. Densities less than 1.003 kg/m? were too difficult 
to measure accurately and with densities above 1.10 kg/m? the solution became opaque 
making visual observation unreliable. 

The compartment was constructed from 6 mm thick acrylic with interior 
dimensions of 0.15 m wide, 0.30 m long (L) and 0.15 m high. Figure 3.2a shows the plan 
and elevation views of the compartment. A flange was built at one end of the 
compartment so that the opening geometry could be easily modified by replacing a face 
plate bolted to the flange. Four opening geometries were used, as seen in Fig. 3.2b, the 
cross hatched area indicates the opening. Opening #1 was a fully open wall, 0.15 m by 
0.15 m, used to demonstrate the similarity between the transient results presented in 
Section 5 and the steady state inviscid theory of Benjamin!. Opening #2 is a horizontal 
slot, 0.15 m wide by 0.05 m high, centered vertically in the end wall and corresponds to 
the opening used in the backdraft experiments. Opening #3 was a 0.05 m square, centered 
vertically and horizontally in the wall designed to simulate a window opening. Opening #4 
was 0.12 m high by 0.05 m wide centered horizontally with the bottom of the opening at 


floor level, to simulate a door. Openings 1 & 2 can be considered two-dimensional in the 
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large scale whereas openings 3 & 4 are clearly three-dimensional. The end opening was 
covered with a vertical sliding partition that was removed to start the experiment. A 
neoprene gasket was placed between the face plate and sliding partition to prevent leakage 
before the experiment began. The compartment was made negatively buoyant by adding 
1.6 kg of lead shot in ballast channels beneath the compartment as indicated in Fig. 3.2a. 
The solution in the compartment was regular tap water with: pH 6.8, density 1.000 
kg/m3, and temperature 18°C. 

In the early two dimensional experiments, blue vegetable dye was added to the 
compartment. A small amount phenolphthalein (4 x 10-°M concentration) was also added 
to the compartment fluid to visualize the gravity current mixing. When phenolphthalein 
mixes with a base, in this case sodium hydroxide, the product of the reaction is red. This 
reaction is believed to be diffusion limited. The red product is strongly visible even in 
dilute concentrations. Turbulence is unaffected by the reaction since it produces little 
surface tension, buoyancy, or heat release. Although reversible, the disappearance of the 
red product can be minimized by keeping the pH in the large tank high, in this case 11.7. 
The pH in the large tank was raised using sodium hydroxide crystals. Unlike the passive 
scalar techniques, such as dye, the chemical reaction of the phenolphthalein gives a much 
better indication of the mixing that occurs in the gravity current. A thorough discussion of 
this technique is given by Breidenthal®. 

In the three dimensional experiments, a mirror was placed above the 
compartment at a 45° angle to show the plan view of the gravity current in the same plane 
as the elevation view for video recordings. When the mirror was used, the blue dye was 
eliminated and the phenolphthalein concentration was micreasediby a factor of 4 to 1.6 x 


10~*M to produce a more visible gravity current. 
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Figure 3.2a - Sketch of the salt water compartment showing the elevation and plan views. 


4] 


Coe NJ 


Ve, 


q 
4 


150 mm 


¢ 
oe. 
cS 


150 mm 
COE OCCCP x 
SSO 
ie 


y 


SS 


ES 


Opening #2 


Open] nla 


SS 


4 


Opening #4 


Opening #3 


Sketch of the four opening geometries for the salt water compartment. 


Figure 3.2b - 
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Once the two solutions were prepared, specific gravity, temperature, and pH were 
recorded. The compartment was then lowered into the tank and the partition on the 
compartment was removed within 120 s to avoid leakage effects. Within 0.1s the partition 
was completely clear of the opening. The gravity current was recorded using a high 8 mm 
video camcorder at 30 frames per second. Typically these experiments lasted less than 


five minutes. 


3.4 Compartment Gravity Current Structure 

Some characteristic features of a steady state gravity current include: a head at the 
front of the current, mixing at the shear interface between the two fluids, and a series of 
advancing lobes and clefts at the leading edge. Figure 3.3 is a simple sketch showing plan 
and elevation views displaying some of these features on a steady state gravity current. 
The foremost point of the current is slightly raised above the bottom surface to a height of 
h,. This lifting of the head is a result of the faster moving heavier fluid overrunning the 
slower light fluid. The pahten fluid is forced under the gravity current as a result of the no 
slip condition at the lower bounding surface. 

The overrun fluid causes a gravitational instability which is largely responsible for 
the three-dimensional effects which are seen in natural gravity currents. The instability is 
manifested as the lobes and clefts which make up the leading edge of the current and the 
billows which form above and behind the head of the current?!°. In Fig. 3.3, the plan view 
of the leading edge shows the lobe and cleft structure. The width of the lobe is b~O(h,) 
as reported by Simpson’. As a lobe widens it will split and a portion of the dense fluid, 
mixed with the lighter fluid overrun by the current, is swept up and over the head forming 
a new billow behind the head of the advancing current. As a result of the split two smaller 
lobes are formed and a cleft develops between them. The billows which form behind the 


head of the gravity current are both qualitatively and quantitatively similar to the Kelvin- 
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Helmholtz instability of a shear layer separating two flowing fluids of different 
densities®!!. 

Figure 3.4a is a photograph showing the gravity current approximately 3L/4 into 
the compartment for the full open condition. The lower half of the photograph shows the 
profile of the gravity current. The top half of the photograph shows the plan view 
reflected in the inclined mirror. The mixed region, confined to a shallow layer between the 
two fluids along the shear interface, is red due to the chemical reaction of the 
phenolphthalein but appears gray in the black and white photograph. In profile the gravity 
current head is raised above the lower boundary as described above. The fluid that is 
overrun by the current is visible in the photograph as the gray area under the head. Along 
the interface between the two fluids the large billows can be seen developing behind the 
head. In plan view the lobes and clefts which make up the leading edge can be clearly 
seen. The lobes and clefts make the leading edge of the current difficult to determine and 
contribute significantly to the overall error analysis described in Section 3.5. 

The photograph in Figure 3.4b of the gravity current approximately 3L/4 into the 
compartment for the slot opening shows an important conclusion from this salt water 
modeling. The current is mixed throughout, as indicated by the increased size of the gray 
(red) region compared to the full open case. The increased mixing is a result of the 
rearward facing step caused by the opening being placed above the floor. Otherwise, the 
profile of the gravity current shows a similar structure to the traditional gravity currents 
discussed above i.e., the slightly raised head and the billows formed behind the head. The 
plan view of the current clearly shows the presence of lobes and clefts at the leading edge. 
The gravity current head is not as high as in the full open case and the lobes are also 
smaller. The retarding effect of the no slip boundary condition along the walls is also 


apparent in the plan view. 
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Figure 3.3 - Sketch of the gravity current with a no slip condition at the lower boundary. 
Both elevation and plan views are shown. Characteristic features of the gravity current 
are shown including the lobes, clefts, and billows. 
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Figures 3.4c & d are photographs showing the gravity currents approximately 3L/4 
of the way into the compartment for the window and door opening conditions 
respectively. By the time the current has reached the 3L/4 point, the effects of the three- 
dimensional opening are reduced and the gravity current is qualitatively similar to the two- 
dimensional slot opening. The effect of the three-dimensional opening can be seen in the 
series of photographs shown in Fig. 3.5. These photos show how the gravity current 
enters the compartment for the window geometry. The current initially spreads radially 
from the opening but by the time the current reaches approximately L/2, the leading edge 
is moving into the compartment in a similar manner to the slot opening. The three 
dimensional opening will increase the amount of entertainment that occurs as the fluid 
cascades over the edges forming the opening. This entrainment is caused by relatively 
large coherent structures somewhat similar to those occurring in the plume as it exits the 


compartment. 
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Figure 3.4a - Photograph of the gravity current approximately 3L/4 into the compartment 


for the fully open condition with B = 0.018. The grid shown on the model is 25 mm 
squares. 


Figure 3.4b - Photograph of the as current approximately “SLA into the compartment 
for a center h,/3 horizontal slot opening with B = 0.024. 
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Figure 3.4c - Photograph of the gravity current approximately 3L/4 into the compartment 
for a central, h,/3 square, window opening with B = 0.032. 
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Figure 3.4d - Photograph of the gravity current approximately 3L/4 into the compartment 
for a door opening of width h,/3 and height 7h,/9 with B = 0.026. 
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Figure 3.5 Series of four photographs of the gravity current modeling showing the 
developing gravity current as it enters the compartment for the window opening B = 
0.032 at times: a 1.4s,b2.2s,c3.3 s,d4.6s 
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3.5 Results of Salt Water Modeling 

The experimental nondimensional velocity Froude Number, as given in Eq. 7, is 
defined in terms of the gravity current velocity, v = L/t,,, where t,, is the time required for 
the leading edge of the gravity current to reach the wall opposite the opening. Times 
reported for these experiments were obtained from frame by frame analysis of the video 
recordings at 30 frames per second. 

Once the gravity current reaches the rear wall it is reflected up and around until it 
travels toward the opening. A Froude number is also calculated for the returning current. 
is the time from 


The returning gravity current velocity is, v, = (2L+2h,/3)/t,,,, where t 


out? out 


opening to the time the reversed current returns to the opening wall. The 2h,/3 factor is 
used to account for the length the current travels up the wall opposite the opening and 
compares well with video observation of the distance the current stretches up the wall. 
Using this factor results in v* being nearly equal for both the entering and exiting currents. 
For the full opening condition it is not possible to determine the leading edge of the 
returning current because there is no restriction on the exiting flow. 

In Table 3.1 8, t,,, v, v*, and Re for the entering wave are given in columns | thru 


5 respectively. For the exiting current, t and v., are given in columns 6 thru 8 


out Veo 
respectively. Looking at the v* values given in Table 3.1 it can be seen that the value is 
constant, differing only with geometry, over the range of B investigated. Figure 3.6 is a 
plot of v* versus B for the entering wave for all four opening geometries. The average 
values of v* and h* = h)/h,, are given in Table 3.2. The fully open value of 0.44 is 
confirmed by independent two dimensional computations presented in chapter 4. This 
value also compares well to the 0.47 < v* < 0.50? reported for lock exchange problems 
and v* =0.5 derived in section 2. The experimental v* = 0.44 is lower than the v* = 0.5 


derived in section 3.2, due to the mixing and the transient flow in this compartment. The 


average v* values for the slot, door and window decrease as the mixing increases. As the 
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slot, door, or window size relative to h, decreases v*—0 from the Table 3.2 values. 
Similarly, increases in relative size cause v*—>0.44. As the aspect ratio (L/h,) increases 
this limit may approach 0.50. 

The error bars shown in Fig. 3.6 were calculated by compounding the errors for 
each parameter in Eqs. 1, 7 and 10. The large relative error bounds for the window 
Opening are a result of the reduction in v*. The absolute error remains unchanged. The 
nondimensional height of the entering gravity current head, h* in Table 3.2, is based on 
visual observation from the video recordings of the experiments. The grid on the 
compartment seen in Figs. 3.4 and 3.5 is used to determine the heights. The grid lines are 
25 mm apart and can be visually divided into four equal parts giving an accuracy of + 
h,/24. 

The average head height of the gravity current is viewed over the distance 3L/4 to L to 
reduce any effects caused by the opening. The average head height did not change over 
the length of 3L/4 to L, within the accuracy of the measurement. The fully open head 
height, h* = 0.5 is consistent with the Eq. 5. The decreasing nondimensional head height 
with increasing mixing is also consistent with the h*<0.5 suggested by Benjamin! for 
dissipative flows. 

The Reynolds number shown in column 5 of Table 3.1 ranged from 


939<Re<13407. The Reynolds number is defined as: 


Re=——-. (3.10) 


In order to compare with other gravity current results, hy is used in this definition. Over 
this range of Reynolds number the nondimensional velocities were found to be constant 
for each opening geometry. This Reynolds number independence is consistent with the 
results of Keulegan!? and Barr!3 who found that the nondimensional velocity was strongly 


dependent on the Reynolds number for Re < O(10%) and independent for large Re. 
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Abraham and Vreugdenhil!* indicate a slight increase in the nondimensional velocity for 
large Re. Simpson and Britter!5 indicate that the nondimensional velocity is either 
independent or only slightly dependent on the Reynolds number for Re>O(103). For a 3m 
compartment fire, which is a candidate for a backdraft, the expected Reynolds number 
range would be 5 x 103<Re< 5x 104. The independence suggested in the literature, and 
shown in Fig 3.6 and Table 3.1 for 103 < Re < 10¢ indicates that the salt water results are 


directly applicable to typical fire compartments. 
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Table 3.1 Summary of the salt water modeling results for the entering and exiting current. 
Results are shown for all four opening: full, slot, door, and window. 
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Figure 3.6 - Nondimensional velocity versus density difference ratio for four opening conditions. 
The O,¢, 0, a represent the data for the full, slot, window and door opening. The (— — -), (-------- 


),(— ---), (- -—-), represent the average values for the full (0.44), slot (0.32), door (0.35) and 
window (0.22), geometries. 


Table 3.2 - Average values for V* and h* calculated from the data reported in Table 3.1. 
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3.6 Conclusions 

This experimental work shows that the gravity current entering a compartment is 
both qualitatively and quantitatively similar to other naturally occurring gravity currents. 
The nondimensional velocity for the full opening compares well with the perfect fluid 
theory presented in section 2. The structure of the entering gravity current head for the 
full open condition shows a shallow mixed region riding on the current at the interface. It 
also has the detailed features of the steady state gravity currents, i.e. billows, lobes, and 
clefts. Entering currents for the slot, window, and door openings show a different gravity 
current structure with the mixed region occupying nearly the entire current due to the 
enhanced mixing near the opening. Similar detailed features appear on the these currents. 

The values of v* and h* obtained here for a variety of opening geometries, are 
independent of the density difference ratio, B. The exiting gravity current is also 
independent of B, and has a v* approximately equal to the entering current. These results 
can be applied to predict the time of ignition for a backdraft with compartment and 
opening geometries similar to the conditions reported here. 

Additional research is necessary to investigate other compartment and opening 
geometries. Future work should focus on using larger aspect ratios (L/h,) and a variety 
of opening ratios (side/h,) as well as openings offset from the wall center line. More 
sophisticated instrumentation may also be used to measure concentrations within the 


current. 
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CHAPTER 4 


NUMERICAL AND EXPERIMENTAL GRAVITY CURRENTS RELATED TO 
BACKDRAFTS 


4.1 Introduction 

When a fire occurs in a closed compartment where the only ventilation is due to 
leakage the fire can become limited by the available oxygen and produce large amounts of 
unburned fuel. Ifthe leakage rate is low enough the fire may enter a smoldering stage. 
Temperatures within the compartment will be low compared with flashover temperatures 
but significantly higher than the ambient. The higher temperatures decrease the average 
compartment density below the ambient density. When the compartment is opened a 
density driven flow referred to as a gravity current enters the compartment. The dense 
ambient air pours through the opening mixing with the hot compartment gas as the current 
travels across the floor. If the fuel concentrations are high enough and the gravity current 
comes in contact with and ignition source, a backdraft will np through the compartment 
injuring any unsuspecting firefighters trapped in the wake of the explosion. A better 
understanding of gravity currents which enter the compartment is necessary to improve 
our knowledge of the backdraft phenomena. 

A large body of research is available on the subject of gravity currents.!2 Much of 
this research involves laboratory scale salt water models. Typically, salt water models of 
gravity currents do not provide complete answers to gravity current questions. Some 
important detail is lost when a salt water model is used including the large vortical 


structure typically seen in shear flow problems. 
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Here, a two-dimensional numerical simulation is used to analyze the gravity 
current as it enters a fire compartment. The results of the two dimensional computations 
provide valuable insight into the detailed structure of the gravity current. The results of 
the numerical simulation are validated by comparison with the results from two different 
experiments. In the first set of experiments, a salt water model with two fluids of different 
density is used to visualize the flow into the compartment as described in chapter 3. The 
second set of experiments uses a large scale, 1.2m x 1.2m x 2.4m, compartment to 
produce backdrafts in the laboratory, see chapter 2. Data collected in the highly 
instrumented backdraft experiments on the entering gravity current are compared with 


both the salt water and numerical results. 


4.2 Numerical Modeling 
Consider the buoyancy-driven flow induced by the interaction of salt water and 
fresh water initially separated by a vertical interface. The equations of motion for this 


incompressible, isothermal mixture are: 


div(u) = 0, (4.1) 
O(pc)/Ot + div(pcu) = div(pDVe), (4.2) 
p(du/ot + uVu) + Vp - pg = pvVu, (4.3) 


where c is the salt water mass fraction (defined as the ratio of the mass of salt water to the 
total mass of fluid in a given volume element) u, is the velocity vector, p is the pressure, g 
is the acceleration of gravity, v is the kinematic viscosity, and D is the diffusion 
coefficient. The latter two quantities will be assumed to be constant. The density of the 


mixture may be expressed as p= p,(1+Pp), where B = (p, —p,)/p,; and p, and p,are the 
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densities of the fresh and salt water, respectively. In terms of B and 0, the mass fraction c 
may be written as: 
c=p(1+B)/(1+Bp). (4.4) 

We are interested here in the motion of the fluid mixture in a two dimensional 
polygonal configuration, consisting of a small chamber initially filled with a fresh water, 
separated from the salt water outside by a vertical interface. Equations (4.1)-(3) are 
solved numerically in nondimensional form using finite differences. The Reynolds and 
Schmidt numbers that result from the scaling are given by Re=vj,/v and Sc=v/D where 
v,is a characteristic velocity related to the Froude scaling v, = VBgh, , and h, is the height 
of the enclosure. An alternating direction implicit numerical scheme was used to solve the 
above equations. Details of the numerical method may be found in Ref. 3. For the 
comparisons with the salt water and backdraft experiments shown here, computations 
were performed on the IBM RISC System/6000 Model 550 of the Mathematics 
Laboratory and the Building and Fire Research Laboratory at NIST. A typical 
computation required between 20 and 80 megabytes of memory and 5 to 25 hours of CPU 
time. The resolution of the computational grid determines the maximum Reynolds number 
for a given run. Roughly, this maximum value scales as K”, where K is the number of grid 
cells in the direction of the length scale h,. The largest Reynolds number reported here is 


50,000, and this simulation required a grid of dimension 1024 x 256. 


4.3 Salt Water Experimental Apparatus and Procedures | 

Salt water experiments were conducted by placing an acrylic chamber within a 
larger glass tank. The chamber was constructed from 6 mm thick acrylic with interior 
dimensions of 0.15 m wide, 0.30 m long and 0.15 m high. Figure 4.1 shows the plan and 
elevation views along with the opening geometries for the chamber. Two two- 


dimensional opening geometries were used, as seen in Fig. 4.1, the cross hatched area 


a 


indicates the opening. The end opening was covered with a vertical sliding partition that 
was removed to start the experiment. The large tank, 0.3 m wide, 0.6 m long, and 0.45 m 
deep, contained a dense water and salt solution ranging in density from 1.003 kg/m? to 
1.101 kg/m3, 

The solution in the chamber was regular tap water with a pH of 6.8 and a density 
of 1.000 kg/m3- A small amount of phenolphthalein (>2x10~M) was added to the 
chamber. Phenolphthalein, a common pH indicator, was used to visualize the gravity 
current. When phenolphthalein mixes with a base, in this case sodium hydroxide crystals 
were added to the large tank to raise the pH to 11.7, the product of the reaction is red. 
This reaction is believed to be diffusion limited. The red product is strongly visible even in 
dilute concentrations. Turbulence is unaffected by the reaction since there is little surface 
tension, buoyancy, or heat release produced by the reaction. Unlike passive scalar 
techniques, such as dye, this chemical 
reaction is a much better indicator of the mixing within the gravity current. A formal 
discussion of this technique is given by Breidenthal?. 

Once the two solutions were prepared, specific gravity, temperature, and pH were 
recorded. The chamber was then lowered into the tank and the partition on the chamber 
was removed within 120 s to avoid leakage effects. Within 0.1s the partition was 
completely clear of the opening. The gravity current was recorded using a high 8 mm 
video camcorder at 30 frames per second. A more complete discussion of the apparatus 


and procedures used in the salt water modeling can be found in chapter 3. 
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Figure 4.1 - Sketch of salt water chamber showing the elevation and plan views as well as 
opening geometries. 


4.4 Backdraft Experimental Apparatus 
A second series of experiments were conducted using a half scale compartment 
filled with hot gases from a methane fueled gas burner. Figure 4.2 is a sketch of the 


compartment. In one of the short walls, a 0.4 m high, 1.1 m wide opening was centered 
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vertically. A 0.3 m square gas burner, with a horizontal surface 0.3 m off the floor and 
centered horizontally along the wall opposite the opening, provided the initial fire within 
the compartment. Gas flow to the burner varied from 70 to 200 kW of technical (98% 
pure) methane. The burner was ignited with an electric arc located 50 mm above the 
burner and centered along the burner edge facing the room (SPK2). The backdrafts were 
ignited using either SPK2 or another spark generator placed 0.15 m above and centered 
over the burner surface (SPK1). The sparks were generated using a 10,000 VDC 
transformer and two 3 mm diameter 308 stainless steel electrodes 5 mm apart. 

A computer controlled hatch, hinged at the bottom, covered the slot until a 
predetermined opening time was reached. In the slot, six bidirectional probes were 
installed in a vertical rake. The probes were 17 mm in diameter and designed in 
accordance with the guidelines given by McCaffrey and Heskestad.> The probes were 
evenly spaced, 65 mm apart centered horizontally in the opening and numbered 
sequentially from the top down. The top (#1) and bottom (#6) probes were 43 mm from 
the soffit and sill, respectively. The bidirectional probes were oriented horizontally to 
avoid incorrect readings due to buoyancy effects. The pressure differential was measured 
using a differential pressure transducer which had a calibrated range of + 25 Pa. The 
response time of the transducer was 30us. Readings were recorded approximately 50 
times a second and a 11 point smoothing routine was applied to the data. Bare bead 
thermocouples made from 0.5 mm type K wire, with an average bead diameter of 1.1 mm 
were placed at each probe to measure the temperature as required for density 
calculations. No correction was applied to the thermocouple data. Aspirated 
thermocouples placed in close proximity to the probes indicated that a radiation correction 
was unnecessary. The response time of these thermocouples was adequate for the gravity 
current, but much too slow to characterize backdraft velocities. Additional details can be 


found in the references 1 and 8. 
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In each experiment, the gas burner was ignited in the closed room. The fire burned 
until the burner could no longer support combustion due to a lack of oxygen. The burner 
was then left on to allow unburned fuel to accumulate. At a predetermined time, the 
burner was turned off and 5 s later the hatch was opened to allow the gravity current to 
enter. In some of the experiments, the spark ignitor at the burner edge (SPK2) was 
activated continuously and when the gravity current reached the spark a backdraft 
occurred. In other experiments, SPK2 was not used and the spark above the burner 
(SPK1) was activated at opening. When the gravity current reached SPK1, ignition of the 
backdraft occurred. A video recording of each experiment was used to determine event 
timing. Video tape and computer data times are synchronized using a computer controlled 


light, visible in the video frame, which turned on when the hatch was activated. 
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Figure 4.2 - Sketch of the half scale backdraft compartment showing important features of 
the apparatus. 


4.5 Qualitative Results 

Figure 4.3a shows the density field from the numerical simulation for the fully open 
condition. The black color represents the compartment fluid and the lightest gray color is 
the ambient fluid. Looking at Fig. 4.3a the gravity current can be divided into two 
regions: region 1 is purely ambient, cold, fluid which is moving toward the head of the 
gravity current and region 2 is the mixed layer along the shear interface which consists of 
hot, fuel rich, compartment fluid and rolled up within long coherent structure of cold, 


oxygen rich ambient fluid. 
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Figure 4.3b is a photograph from the salt water experiments for the fully open case 
at approximately the same location as Fig. 4.3a. The photograph closely resembles the 
numerical simulation results. Region 1 1s the clear salt water which indicates that no 
phenolphthalein has reached that portion of the gravity current. The mixed layer, region 2, 
is the gray (red) area along the shear interface. The large vortical structure seen in the 
numerical simulation is not visible here because the photograph of the salt water 
experiment is an integral along a line of sight across the entire width of the chamber and 
the details of the structure are lost. Experiments by Simpson® show similar large vortices 
when vertical slit lighting and fluorescent dye are used to illuminate a more two 
dimensional image of a steady state gravity current. 

Figure 4.4a shows the computed density profile for the h,/3 centered slot opening 
condition. The structure of the gravity current is significantly different from the full 
opening case. The current can not be divided up into two distinct regions as in fully open 
case. Large vortices make up entire gravity current for the slot opening. The coherent 
structure making up the gravity current indicates that large scale mixing is occurring. The 
increased mixing is caused by the h,/3 centered slot opening acting as a rearward facing 
step which is a well known source of large vortical structures. 

Figure 4.4b is a photograph of the salt water experiment with the h,/3 centered slot 
opening showing the gravity current in approximately the same location as Fig. 4.4a. The 
large 
scale mixing predicted by the computations is seen as the dark gray color throughout the 
gravity current. From these results, it can seen that the simple two region model used for 
steady state gravity currents works for the full opening but cannot be applied to the slot 


opening condition. 
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Figure 4-3a - Density profile for the full opening case. Shown here after 4.0 
nondimensional time units. Reynolds number is 20,000. 
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Figure 4.3b - Photograph of the gravity current approximately 3L/4 into the compartment 
for the fully open condition. The grid shown on the model is 25 mm squares, B= 0.018. 
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Figure 4.4a - Density profile for the h,/3 centered slot opening case. Shown here after 6.8 
nondimensional time units. Reynolds number is 20,000 
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Figure 4.4b - Photograph of the gravity current approximately 3L/4 into the compartment 
for h,/3 centered slot opening, B= 0.024. 
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4.6 Quantitative Results 

To quantitatively compare the numerical simulation with the experiments, the 
gravity current transit time is used. The transit time is the time required for the leading 
edge of the gravity current to reach the wall opposite the opening. For the numerical 
simulation the transit time was determined when the density at the rear wall changed by 
10%. In the salt water experiments, the transit time was taken from video recordings of 
the gravity current. In the half scale backdraft experiments, the gravity current is not 
directly measured. The time to ignition is used to approximate the transit time assuming 
that the gravity current is ignited as soon as the current reaches an ignition source. 

Figures 4.5a and b show the transit times versus relative-density difference, B, for 
the full and slot opening conditions, respectively. Figure 4.5a shows the excellent 
agreement between the salt water model and the numerical simulation for the fully open 
condition. For the h,/3 centered slot opening, Fig. 4.5b, the agreement is also good but 
computed values are consistently longer than the values measured in the salt water 
experiments. 

Figure 4.6 shows the comparison between the time to ignition from the half scale 
experiments and the numerical results versus 8B. The o indicates that ignition occurred at 
the SPK2 and the o indicates that ignition occurred at the spark above the burner, SPK1. 
Although the data are somewhat scattered, the comparison indicates the expected trend. 
For the burner spark ignition the results are excellent. The longer times seen for the spark 
above the burner maybe due to the burner's effect on the flow field. The two data points 
which occurred earlier and the three which occurred later than expected may demonstrate 
the dependence of the ignition phenomenon on the steep concentration gradients shown in 


Fig. 4.4a. 
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Figure 4.5a - Plot of the transit time versus density difference, B, for the full open 
condition comparing the numerical simulation ( -) with salt water modeling results (0). 
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Figure 4.5b - Plot of the transit time versus density difference, B, for the h,/3 centered 


slot opening condition. Compares the numerical simulation (-----) with salt water 
modeling results (0). 
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Figure 4.6 - Plot of the time to ignition for the half scale backdraft experiments and 
numerical simulation transit time versus density difference, B. The (0) indicates ignition at 
the burner spark, (0) indicates ignition at the spark above the burner, and (. ) indicates 
the computational results. 


A more direct comparison between the numerical simulation and the half scale 
backdraft experiments can be made using the bidirectional probe measurements. The 
probe velocity is calculated using the relationship suggested by McCaffrey and 
Heskestad:? 

EIGER) = (4.5) 


where v is the velocity at the probe, Ap is the pressure difference measured across the 
probe, p is the local density, and C(Re) is an empirical calibration constant which is a 
function of the Reynolds number. Typically, C(Re) is taken as a constant of 0.926 which 


gives a maximum error of 7% for Re > 520°. In this case, the Re was smaller, Re ~ 300, 
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so the following calibration constant was required: 


po em 
C(Re) 


= 1.533-1.366x 107 Re+ 1.688 x 10° Re? —9.706x 10" Re? +2.55x 107? Re*—2.484 x 107 Re® 
Equations (4.5 and 4.6) are applicable for 40 < Re < 3800. For Re < 40, v<0.07 m/s 
which is considered negligible. 

Figures 4.7a and b compare the velocity history for probes 1, 2, 5, & 6 with the 
numerical simulation for B=0.52. Negative velocities indicate flows into the compartment. 
The hatch was released at 0 s. At approximately 1.5 s the opening hatch strikes the table 
and causes excessive noise in the data. From 0 to 4 s the flow in the opening is 
developing. After approximately 4 s the flow can be assumed to be quasisteady. At 4 s 
the total mass flow into the chamber is 0.28 kg and the gravity current is approximately 
3L/4 into the compartment as shown in Fig. 4.3. If ignition does not occur, the flow will 
slowly diminish to zero as the compartment cools. The flow into and out of the | 
compartment which continues long after the initial gravity current has subsided is the 
result of the thermal energy stored in the chamber heating the incoming air and driving the 
flow. In one experiment where the ignition was delayed for 600s, the velocity dropped to 
~0.2 m/s. 

Figure 4.8 shows the vertical velocity profile in the opening at different times. The 
lines show the numerical simulation and the symbols are the experimental data. At 4s the 
data indicate that the flow has developed into the expected in/out profile and show 
excellent agreement with the numerical simulation. The velocity profile shows little 


change at 8 s when the spark is activated and ignition occurs. 
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Velocity (m/s) 


Time (s) 


Figure 4.7a - Plot showing the velocity history for the probe #1 (0) and #6 (o)in the 
opening of the half scale backdraft compartment compared with the numerical simulation 


results shown as, probe #1 (-----) and probe #6 ( ). 


Velocity (m/s) 


Time (s) 


Figure 4.7b - Plot showing the velocity history for the probe 2 (0) and 5 (0) in the 
opening of the half scale backdraft compartment compared with the numerical simulation 


results shown as, probe 2 (-----) and probe 5 ( ). 
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Figure 4.8 - Plot of the velocity profiles in the opening of the backdraft compartment 
shown at 4 s and 8 s after opening. Numerical simulation results are shown as lines and 
the experimental data is shown as symbols. (4 s O-----) (8s 0 ) 


Figures 4.9a and b are plots of quasisteady velocity in the opening versus B. Data 
are shown as discrete points, numerical simulation results are the solid lines, and potential 
flow results are the dashed lines’. There is good agreement between the experimental data 
and the numerical results, with the numerical result being consistently higher than the data. 
The potential flow results is almost a factor of 2 greater than the measured result which 
indicates that assumption, as expected, oversimplifies the problem. The numerical 
simulation is much closer to the measured values. The differences between experiments 
and computation are likely to be due to turbulence effects not included in the simulation 
and the large errors which can be expected in the measurements at this low velocity. For 


velocities of O(1 m/s), the experimental error is as high as +40%. 
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Figure 4.9a & b - Plot of the quasisteady velocity in the opening of the backdraft 
compartment versus density difference ratio, B. The experimental data are indicated as 
symbols, the numerical simulation is shown as solid lines, and the potential flow result is 
shown as a dashed line. (probe 1 o)(probe 2 0)(probe 5 A)(probe 6 +) 
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4.7 Conclusions 

The two dimensional numerical simulation presented here accurately predicts the 
compartment gravity current for the fully open condition. The two dimensional 
computational density profile shows that the main body of the current is made up ambient 
fluid while a mixed layer made up of large vortical structures exists along the shear 
interface between the two fluids. Such structure is also observed in the salt water 
experiments although the individual vortices are not resolved due to the limitations of the 
flow visualization techniques. Results of other researchers!° indicate that large vortices do 
occur behind the head of a steadily propagating gravity current which are similar in 
appearance to the transient compartment gravity current results shown here. 

For the slot opening, the results of the two dimensional simulation still compare 
very favorably with the experimental results. The computed density profile shows that the 
structure of the gravity current for the h,/3 centered slot opening is significantly different 
from the fully open condition. For the slot opening, the entire gravity current is filled with 
a complex structure of large scale vortices which translates into large scale mixing. This 
large scale mixing for the h,/3 centered slot opening results from the rearward facing step, 
formed by the lower edge of the slot. 

Transit times predicted by the numerical simulation compare well with the salt 
water experiments although the computed times are slightly longer for the slot opening. 
Comparing the computational transit time with the time to ignition for the backdraft 
experiments gravity current also shows good agreement. Numerical velocity profile 
predictions in the opening compared well with backdraft experiment results, falling well 
within the experimental error bounds. Future work should focus on different 
compartment aspect ratios, smaller openings, and openings at the floor level. Improved 


flow visualization techniques including slit lighting to obtain a two dimensional image of 


75 


the gravity current are recommended. Time dependent concentrations should be measured 


within the gravity current for comparison with computational results. 
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CHAPTER 5 


QUANTITATIVE BACKDRAFT EXPERIMENTS 


5.1 Introduction 

The dangerous consequences of a backdraft are documented in numerous fire 
service publications and training manuals.!23-4 However, little research has been done in 
the area of backdrafts, and only recently has a scenario been presented describing the 
fundamental physics underlying backdraft phenomena. A backdraft is defined as a rapid 
deflagration following the introduction of oxygen into a compartment filled with 
accumulated unburned fuel. 

The scenario presented here assumes a fire ina closed room. The fire heats up the 
room, and leakage in the bounding surfaces minimize the pressure differential. The hot 
layer descends over the fire as the oxygen concentration is reduced and the combustion 
efficiency decreases. Excess pyrolyzates accumulate in the upper layer forming a fuel rich 
mixture of low oxygen content. A small flame or glowing ember exists as a source of 
ignition. Suddenly, a new ventilation opening is provided and cold, oxygen rich, air enters 
the compartment and propagates across the floor as a gravity current. Large scale mixing 
in the gravity current provides areas within the flammable range which can ignite when 
they contact a source of ignition. Once ignited, a flame propagates through the 
compartment and drives the remaining unburned fuel out through the opening to burn 
outside the compartment in a spectacular fireball. 

In this paper, experimental results are presented from a series of half scale 
experiments attempting to quantify backdraft. Experimental variables included fuel flow 


rate, burn time, ignition location, ignition delay time, burner height, species sample 
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location, and opening size. This paper focuses on 17 experiments in which the opening 
geometry was a vertically centered, horizontal slot, in one wall; the fuel source was a 0.30 
m square burner 0.30 m above the floor; and the ignition source was a spark located 
opposite the opening. Two different burner flow rates were used, 70 kW and 200 kW. 
Data collected in these experiments had two goals: 1) to characterize the conditions in the 


compartment prior to backdraft and 2) to quantify the severity of the deflagration. 


5.2 Experimental Design & Procedures 
5.2.1 Apparatus: Experiments were conducted in a special compartment designed to 
safely control the dangerous overpressures expected in backdrafts. The experimental 
apparatus dimensions were limited to half a small residential room to minimize the 
expected hazard and to allow the experiments to be conducted inside a 900 m? facility. 
Figure 5.1 shows a schematic of the apparatus giving the internal dimensions of the 
compartment and the locations of the instrumentation. In order to control the 
overpressure hazard, one long wall was a pressure relief panel designed according to refs. 
5 and 6. The interior surfaces of the compartment were lined with a 50 mm thick 
refractory fiber blanket installed over the gypsum wallboard on the walls and ceiling to 
provide the primary thermal resistance for the structure. This insulation allowed for 
repeated experiments without the need to rebuild the compartment. A 0.9 m high by 1.5 
m wide observation window of Neoceram’ was installed in the wall opposite the pressure 
relief panel. 

To simulate a window or door, a 0.4 m high by 1.1 m wide opening was centered 
in the short wall opposite the burner, see Fig. 5.1. This opening was covered with a 
computer activated hatch which was opened after the fire had been burning for several 
minutes. A methane burner, 0.3 m square and 0.3 m high, was used in all of these 


experiments. The burner was placed against the wall opposite the opening, as seen in Fig. 
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5.1. A pilot flame was used to ignite the burner and was turned off 10 s after the start of 
the experiment. The primary ignition source for the backdraft was a spark ignitor located 
0.45 m above the floor and centered over the top of the burner. For 3 of the experiments, 
the primary spark malfunctioned and the backup spark ignitor used to ignite the pilot light 
had to be used to ignite the backdraft. The backup spark was located 0.35 m above the 
floor and centered on the side of the burner facing the opening. A 10,000 volt transformer 
was used for each spark ignitor to produce an arc between two 3 mm diameter 308 


stainless steel electrodes 5 mm apart. 
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Figure 5.1 - Sketch of the half scale backdraft compartment showing important features of 
the apparatus. 
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Every effort was made to seal all construction holes to control leakage. The 
primary source of leakage into the compartment was found to be around the pressure 
relief panel and the opening hatch. Gaskets made from the refractory fiber blanket were 
compressed around the edges of these opening to reduce the leakage. A small 0.1 m 
diameter pressure relief vent was placed at the floor level to relieve the pressure from the 
initial burner ignition. Without this vent, a pressure rise sufficient to activate the pressure 
relief panel was produced. A computer controlled cover closed over this vent 15 s after 
ignition. Additional details of the apparatus can be found in Appendix A. 

5.2.2 Species Concentration: In order to characterize the compartment conditions prior 
to a backdraft, the species concentration histories in the upper layer were recorded. Gas 
concentrations measured were: oxygen (O,), carbon dioxide (CO,), carbon monoxide 
(CO), and total hydrocarbons (HC). Continuous gas samples were taken with stainless 
steel probes located 0.6 m from the opening wall, 0.2 m from the ceiling, and 0.6 m from 
the side wall, as shown in Fig. 1. The O,, CO, and CO, samples were taken through an 
unheated sample line in which the soot and water vapor were removed by glass fiber and 
desiccant filters, respectively. 

The HC concentration was more difficult to measure. The hydrocarbon sample 
required a separate heated sample line to prevent the loss of hydrocarbons due to 
condensation. A flame ionization detector was used to measure the hydrocarbon 
(methane) concentration. The effective range for this meter was 0 to 1%. The expected 
range of hydrocarbons was of the order of 20%, by mass. It was therefore, necessary to 
dilute the sample. The dilution system was designed to mix the compartment sample with 
heated ambient air in a ratio of 20 to 1. The flow rate of the dilution air and the sample 
were determined by measuring the pressure drop over a fixed length of tubing. Sample 
flow rates were monitored continuously during the experiment. Typical dilution ratios 


would vary from 22 to 25 depending on the compartment gas temperature. A hot water 
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jacket around the hydrocarbon sample line and dilution air line kept the gases over 60°C, 
well above the maximum calculated dew point of 44°C. To obtain the final species 
concentrations in the upper layer it was necessary to calculate the concentration of H,O 
that was in the upper layer since there was no direct measurement of the H,O. Detailed 
species balances were performed and can be found in Appendix B. 
5.2.3 Temperatures: A vertical thermocouple tree was placed 0.6 m from the opening 
wall and 0.2 m from the pressure relief panel, as shown in Fig. 5.1. The thermocouples 
were made from 0.5 mm type K thermocouple wire with a stainless steel overbraid. The 
average bead diameter was 1.1 mm. The ten thermocouples were located at 0.10 m 
intervals, with the highest thermocouple at 0.15 m below the ceiling. The temperatures 
reported here are uncorrected values. 

The thermal interface height history was calculated from the time dependent 
temperature profiles recorded from the thermocouple tree. The profiles were converted 
into unsteady average upper and lower layer temperatures using the method Quintiere et. 


al.8 applied to steady state temperature profiles: 


hy 
{($oe= [mmr +n, /t, (5.1) 
0 
h, 
JTdx =[h, -h, JT +h, T”, (5.2) 
0 


where T”” and T™ are the upper and lower layer temperatures, and h, and h, are the 
heights of the compartment and the layer interface, respectively. Equation (5.1) is a mass 
balance and Eq. (5.2) retains the same mean temperature as in the data. Assuming that the 


lower layer temperature was the arithmetic average of the two lowest thermocouples 
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allowed the upper layer temperature and thermal interface location to be calculated from 
Eqs. (5.1 and 5.2). 

5.2.4 Compartment Pressure: The compartment pressure history was recorded using an 
electronic pressure transducer with an effective range of 0 to 1250 Pa. The pressure port 
was mounted in the stationary wall opposite the pressure relief panel at floor level. The 
ambient pressure reference was taken outside the building. 

5.2.5 Hatch Flow: The flow in and out of the compartment after the hatch was opened 
was recorded using six bidirectional probes in the h/3 centered slot. The probes were 
located in the horizontal center of the opening and 65 mm apart. The outer probes were 
43 mm from the soffit and sill. Probe velocities were calculated using the relationship 


given by McCaffrey and Heskestad:?: 


v =C(Re) oa (5.3) 


where v is the velocity at the probe, Ap is the pressure difference measured across the 
probe, p is the density at the probe, and C(Re) is an empirical calibration constant which is 
a function of the Reynolds number. Typically, C(Re) can be taken as a constant at 0.926, 
which gives a maximum error of 7% for Re > 520. In these experiments Re < 520 was 
recorded and a Re correction was applied!2. 

The temperature was recorded at each bidirectional probe in the opening using a 
bare bead 24 gauge Type K thermocouple with high temperature glass insulation. Four 
aspirated thermocouples, designed according to Newman et. al.,!° were installed to correct 
the temperatures for radiation. Data from the aspirated thermocouples revealed that no 


radiation correction was required. 


$3 


The mass flow was calculated by integrating the velocity and density profiles over 


the height of the opening as: 


mh, =C f*pvbdy, (5.4) 
where m, is mass flow in the lower part of vent, C 1s the flow coefficient taken as 0.68, b 
is width of opening, and h, and h, are the height of the sill, and neutral axis, respectively. 
5.2.6 Data Acquisition System: Data from each sensor was recorded using a HP 
VECTRA 80486-33 computer with a 8 channel multifunction analog and digital 
input/output board. Two 32 channel analog input muliplexors were connected to this 
system. A total of 31 thermocouple and 17 voltage channels were used. To increase the 
scan rate, data was written to a RAM drive and then down loaded to a file on the hard 
disk immediately following the experiment. The system was capable of recording each 
channel 50 times a second. For experiments greater than 600 s the data was collected at a 
rate of 10 scan/s until 20s before opening when the rate was automatically increased to 50 
scan/s. The reduced scan rate for the initial period of the experiment was done to reduce 
the data file size. Files for these experiments ranged in size from 4 to 10 Mbytes. 

In addition to recording the data, the computer also controlled the experimental 
procedures using solid state relays activated by a digital input/output board also installed 
on the computer bus. The computer controlled systems included, experimental clock, 
burner pilot light, fuel flow, vent cover, spark ignitors, opening hatch, and still camera. 
Each system had a manual override. 

5.2.7 Procedure: Before each experiment a 60 s baseline was taken to record the initial 
conditions. A pilot flame was ignited at the burner 5 s before the start of the experiment. 
At 0 s a solenoid was opened on the methane flow to the burner and the clock was reset to 
zero. The burner was left on for a predetermined time period. Gas flow to the burner was 


terminated 5 s before the hatch was opened. In some of the experiments, a spark was left 
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on the entire time, while in other experiments the spark was not turned on until the hatch 
was opened. When the spark was left on, a dancing flame would often appear and 
consume some of the available hydrocarbons making consistent results difficult to attain. 
In later experiments the spark remained off until the hatch was opened. The dancing flame 


phenomenon is discussed in chapter 2. 


5.3 Experimental Results 

5.3.1 70 KW Fuel Flow Rate: The experimental parameters for the results presented in 
Figure 5.2 are: gas flow rate of 70 kW, burner flow time of 775 s, hatch opening at 780 s, 
and backdraft ignition above the burner. Figure 5.2a shows the upper layer species mass 
fraction histories for O,, CO, CO,, and HC. Idealized layer temperatures and height 
calculated from Eqs. (5.1 and 5.2) are shown in Fig. 5.2b. For the first 80 s, the 
temperature in the compartment rises as the layer within the compartment descends over 
the fire. The O, concentration is dropping as the fire consumes the available O, and the 
CO, concentration increases as a result of the combustion. Over the time period of 80 s to 
170 s the temperature in the compartment drops as the fire is diminished. Over this same 
period the HC concentration is starting to rise as the O, and CO, level off. Waviness seen 
in the upper layer temperature history over the range of 140s to 170 s are the result of the 
flames pulsing before extinction. After 170 s the fire is completely out and the 
compartment begins to cool as seen in the exponential decay in the temperature profile. 
After 240 s the hydrocarbon concentration is steadily increasing. The oxygen 
concentration is slightly increasing as air leaks into the compartment and the CO, 
concentration declines as compartment gases are lost by leakage. At 720 s the HC 
analyzer became saturated and the slope of HC concentration approaches zero. At 780 s 
the hatch is opened and a the gravity current enters the compartment. Once the gravity 


current reaches the ignition source, a flame travels through the mixed region, stirs the 
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compartment, drives combustible gases out the hatch, and culminates in a large external 
fireball approximately 4 m in diameter. The spikes in the temperature shown in Fig. 5.2b 
are caused by the wave propagation through the compartment. Gas concentrations after 
780 s are unreliable due to the highly transient effects of the backdraft. 

The idealized two zone approximation is compared with the temperature data from 
the thermocouple tree at 80 s and 780 s in Fig. 5.2c. At 80 s the compartment 
temperature is at its maximum and there is a substantial temperature gradient in the upper 
layer. The layer is located just above the burner surface. At opening, 780s, the layer is 
still near the top of the burner and the compartment temperature is nearly uniform 
vertically indicating that a one zone approximation is reasonable at this time. Figure 5.2d 
is a photograph from this experiment taken 1.5 s after ignition. Notice the flame burning 


along the top of the gravity current. 
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Figure 5.2a - Species concentration histories from backdraft experiment for O, ), 
CO(— -— -), CO, (------), and HC (——-) for the 70 kW (12th row in Table 5.1) fire source. 
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Figure 5.2b - Idealized two zone upper layer ( ‘), lower layer (— ——), and layer height 
(.-----) histories from a 70 kW (12th row in Table 5.1) backdraft experiment. 


Height (m) 
a) 


300 350 400 450 500 550 600 650 700 750 800 
Temperature (K) 


Figure 5.2c - Idealized two zone approximation compared with temperature data at 80 s 
(-----,0) and 780 s ( ,0). 
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Figure 5.2d - Photograph showing flame structure for the same 70 kW fire source 
backdraft experiment. Photograph taken 1.5 s after ignition of the backdraft. 


Figure 5.2e shows the effects of a backdraft, i.e., short histories for the 
compartment pressure and the total mass flow into through opening. The time starts at 
opening, 780 s and shows a 20 s period during which the gravity current enters the 
compartment, a backdraft is ignited and a large fire ball exits the compartment and 
conditions relax to quasisteady equilibrium. At 780 s the compartment pressure drops as 
the hatch falls open. At ~ 781.5 s the hatch strikes the table and causes a tremor in the 
pressure data. Ignition occurs at 785.6 s, marked with a 1 in Fig. 5.2d. Flames exit the 
compartment at 787.5 s as the pressure reaches the first peak, marked with a 2 in Fig 5.2e. 
The second peak is a result of the large fireball which exits the compartment. Even with 
the large indoor facility and approximately 11 m? of vent area, the pressure rises 
significantly as a large fireball erupts outside the compartment. The large drop after the 
spike is a repercussion of the large fireball. The total mass flow into the compartment is 


calculated from Eq. 4. Negative values indicate flow into the compartment. After 
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opening and prior to ignition, the total mass which flows into the compartment shows a 
steady increase. After ignition, the flow into the compartment is reversed and all of the 
flow is out of the compartment as indicated by the period of zero slope in Fig. 5.2e. Once 
the flame has left the compartment and the fireball has subsided, the mass flow rate into 


the compartment reaches a quasisteady state as indicated by the constant slope in the last 


10 s of Fig. 5.2e. 
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Figure 5.2e - Compartment pressure ( ) and total mass inflow (------) histories for the 
same 70 kW fire source backdraft. Arrow 1 indicates ignition of the backdraft and arrow 


2 indicates flame out the opening. 
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5.3.2 200 kW Fuel Flow: Typical results from a 200 kW experiment are shown in Fig. 
5.3 a-e with the following experimental parameters: burn time of 175 s, hatch opened at 
180 s, and backdraft ignition at the spark above the burner by a continuous spark. Figure 
5.3a shows the species mass fraction histories for the O,, CO, CO,, and HC in the upper 
layer. Idealized layer temperatures and height calculated from Eqs. (5.1 and 5.2) are 
shown in Fig. 5.3b. As the temperature rises in the first 30 s, the O, concentration drops 
and the CO, concentration increases. After 30 s the burning is reduced as the O, reaches 
a minimum, the affects of the reduced burning can be seen in the declining CO, 
concentration and temperature. The HC concentration starts to increase after 30 s and 
increases to approximately 12% at 113 s when a dancing flame ignites at the spark and 
moves around the floor of the compartment consuming hydrocarbons in the upper layer. 
The temperatures in the upper layer rise due to the dancing flame. The slow response time 
of the thermocouple tree causes the temperature rise to occur at 120 s instead of 113 s 
when the dancing starts. The dancing flame stops at 130 s when the oxygen in the lower 
layer is consumed and the hydrocarbons start to build again. 

At 180 s the hatch is opened and a the gravity current enters the compartment. At 
ignition the flame appears to propagate more through the main body of the gravity current 
rather than along the shear interface as seen in the 70 kW case. Then the flame exits the 
compartment and a large fireball approximately 2 m in diameter, considerably smaller than 
the 70 kW fireball, is produced. The spikes in the temperature shown in Fig. 5.3b are 
caused by the flame propagation through the compartment. Gas concentrations after 180 
s are unreliable due to the highly transient backdraft. 

The idealized two zone approximation is compared with the temperature data from 
the thermocouple tree at 25 s and 180 s in Fig. 5.3c. At 25 s the compartment 
temperature is at its maximum and there is a substantial temperature gradient in the upper 


layer. The layer is located just above the burner surface. At opening, 180s, the layer is 
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still near the top of the burner and the compartment temperature gradient is considerably 


reduced. Figure 5.3d is a photograph showing the flame structure 1.5 s after ignition. 
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Figure 5.3a - Species concentration histories from backdraft experiment for O, ( ), CO 
(——-), CO, (------), and HC (—---) for the 200 kW (15th row in Table 5.1) fire source. 
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Figure 5.3b - Idealized two zone upper layer ( -), lower layer (— — —), and layer height 
(------) histories from a 200 kW (15th row in Table 5.1) backdraft experiment. 
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Figure 5.3c - Idealized two zone approximation compared with temperature data at 25 s 
(-----,0) and 180 s ( ,©) for the same 200 kW (15th row in Table 5.1) backdraft 


experiment. 


Figure 5.3d - Photograph showing flame structure for the 200 kW (15th row in Table 5.1) 
fire source backdraft experiment. Photograph taken 1.5 s after ignition of the backdraft. 
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Figure 5.3e shows the effects of a backdraft for the 200 kW fire source, 1.e., 
histories for the compartment pressure and the total mass that has flowed in through the 
hatch since opening. The time starts at opening, 180 s and shows a 20 s period in which 
the backdraft is ignited and the fireball exits the compartment and conditions return to a 
quasisteady. At 180 s the compartment pressure drops as the hatch falls open. At ~181.5 
s the hatch strikes the table and causes a tremor in the data. Ignition occurs at 186.7 s 
marked with a 1 in Fig. 5.3e. Flames exit the compartment at 189.2 s after the pressure 
has reached the first peak marked with a 2 in Fig 5.3e. As seen in the 70 kW case, the 


second peak in Fig. 5.3e is a result of the fireball which exits the compartment. 
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Figure 5.3e - Compartment pressure ( ) and total mass inflow (------) histories for the 
200 kW Figure 5.3d - Photograph showing flame structure for the 200 kW (1 Sth row in 
Table 5.1) fire source backdraft experiment. Arrow 1 indicates ignition of the backdraft 
and arrow 2 indicates flame out the opening. 


5.3.3 Summary: Table 5.1 is a summary of the 17 experiments reported here. columns 1 
and 2 are the burner characteristics, i.e., the burner flow rate and the time the burner gas is 
flowing. Columns 3 - 6 are the compartment species concentrations at opening for O,, 


CO, CO,, and HC, respectively. Columns 7 -9 are results calculated from Eqs. (5.1 and 
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5.2) for the lower layer temperature, the upper layer temperature, and the layer height at 
opening. Column 10 is the observed ignition delay time. Values shown in parentheses 
were ignited at the burner spark, all other backdrafts were ignited at the spark above the 
burner. Column 11 is the peak pressure measured inside the compartment. Column 12 is 
the total mass carried into the compartment after hatch opening. Column 13 is an estimate 
of the diameter of the fireball which burns outside the compartment. Column 14 is the run 
name which corresponds to the data given in appendix C. 

In order to obtain a quiescent environment within the compartment, the gas flow 
rate was set at the low value of 70 kW and the spark ignitors were left off until the hatch 
was opened. The flow rate was more difficult to control at the 70 kW rate and some 
minor fluctuation is seen in the data reported in column 1. Burn times ranged from 295 s 
to 775 s. Times greater than 775 s were felt to be too hazardous to attempt safely. The 
species concentrations are nearly constant for O,, CO, and CO, indicating the repeatability 
of the system. The HC concentrations histories are similar to Fig. 5.2a with a long gradual 
build up, although the slope of the curve changed slightly depending on the burner flow 
rate. The idealized layer temperatures decreased as the burner times increased due to the 
energy loss to the boundaries. Layer height maybe taken as constant throughout the 70 
kW experiments considering the calculation method used. The maximum pressure and 
size of the fireball can be considered as gross measurement of the intensity of the 
backdraft and are seen to increase with the HC concentration as expected. 

_ Experiments were also conducted using a 200 kW fire source and varying the burn 
times from 115 s to 235 s. In all five experiments, the spark ignitor was left on 
throughout the experiment and a dancing flame was observed in all but one experiments. 
The dancing flames causes large thermal instabilities within the compartment and increased 
mixing between the upper and lower layer. The increased turbulence makes it difficult to 


obtain repeatable conditions. 
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5.4 Conclusions 

The results presented here show that the HC concentration must be >10% in order 
for a backdraft to occur. When the HC concentration is < 10% the flame travel is slow 
and the compartment overpressure is much lower. As the HC concentration increases the 
compartment overpressure increases and the backdraft becomes more severe. The 70 kW 
burner flow rate experiments produced excellent backdrafts. The species concentrations 
show long slow changes. The results presented here suggest that the 200 kW source was 
too large for the compartment. Experiments with this large source were highly transient 
and compartment conditions are too unstable to interpret trends. The 200 kW experiments 
were further complicated by dancing flames due to the ignitors being on throughout the 
experiment. 

Future work should concentrate on designing and building a full scale apparatus 
using a controlled fire source and possibly more realistic fuels. Improved gas analysis 
would also be useful to determine what hydrocarbons are present in the real fuel 
experiments. Actual opening geometries could also be used to investigate the effects of 
typical doors and windows. Openings in the ceiling could also be incorporated to study 


the effects of firefighter ventilation tactics. 
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CHAPTER 6 
BACKDRAFT EXPERIMENTS USING A SIMULATED WINDOW OPENING 


6.1 Introduction 

A backdraft is defined as a rapid deflagration following the introduction of 
oxygen into a comnarinen, filled with accumulated unburned fuels. The dangerous 
consequences of backdraft to both firefighters and civilians are well documented in 
numerous fire service publications and training manuals.!34 However, little research has 
been done on the fundamental phenomena underlying backdrafts. 

Consider a fire in a closed compartment where the only ventilation provided is by 
leakage. As the fire heats the compartment, leaks in the compartment bounding surfaces 
permit outflows that minimize any pressure differential’. As the hot layer descends over 
the fire, the available oxygen is reduced and the combustion efficiency decreases. Excess 
pyrolyzates® accumulate upper layer forming a fuel rich upper layer. Suddenly, a new 
ventilation source is provided by a window breaking or door opening. Cold, oxygen rich 
air enters the compartment and propagates across the floor as a gravity current. If an 
ignition source is not immediately available, the gravity current will reflect off the rear wall 
and propagate back to the opening. Once the gravity current has reached the opening, the 
compartment becomes reservoir in a reservoir filling problem. A new lower layer made up 
of a mixture of fuel rich upper layer and oxygen rich incoming air continues to grow. 

Once an ignition source is available, e.g. smoldering ember can be or small flame, the 
backdraft ignites and a deflagration wave propagates through the compartment driving the 
unburned fuel out the opening where it combusts in the spectacular fireball commonly 


associated with backdraft. 
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Experimental results are presented from a series of 11 experiments conducted in 
the half room scale apparatus shown in Fig. 6.1. Here, where the opening geometry is a 
0.4 m square window centered, on a short wall, a dramatic horizontal flame jet exits the 
opening prior to the formation of the fireball. Experimental variables include fuel flow 
rate, burn time, ignition location, ignition delay time, burner height, species sample 
location, and opening size. The fuel source was methane from a 0.3 m square sand 
burner centered at 0.3 m height along the wall opposite the opening. The ignition source 
for the backdraft was two sparks located 0.15 m in front of the burner. Data are 
presented which characterize the conditions in the compartment prior to backdraft and 


quantify the severity of the deflagration. 


6.2 Experimental Design & Procedures 

6.2.1 Apparatus: Experiments were conducted in a special compartment 
designed to safely control the dangerous overpressures expected in backdrafts. The 
experimental apparatus amit were limited to half a small residential room to 
minimize the expected hazards. It also allowed the experiments to be conducted inside a 
900 m3 building at the Richmond Field Station of the University of California at Berkeley. 
Figure 6.1 shows a schematic of the apparatus giving the internal dimensions of the 
compartment and the instrumentation locations. In order to control the overpressure 
hazard, one long wall was a pressure relief panel’®. A 0.9 m high by 1.5 m wide 
observation window of Neoceram? was installed in the wall opposite the pressure relief 
panel. 

To simulate a window breaking due to thermal stress, a 0.4 m high by 0.4 m wide 
opening was centered in the short wall opposite the burner, as shown in Fig. 6.1. This 
Opening was covered with a computer activated hatch which was opened after the fire had 


been burning for several minutes. The burner was placed against the wall opposite the 
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opening, as seen in Fig. 6.1. The ignition source for the backdraft was two spark ignitors 
0.15 m in front of the burner and 0.15 m and 0.3 m above the floor. A 10,000 volt 
transformer was used for each spark ignitor to produce an arc between two 3 mm 


diameter 308 stainless steel electrodes spaced 5 mm apart. 


Pressure Relief 
Panel 


Opening Hatch 


Ze 


Figure 6.1 -Sketch of the half-scale backdraft compartment showing important features of 
the apparatus. 
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Every effort was made to seal all construction holes to control leakage. The 
primary source of leakage into the compartment was around the pressure relief panel and 
the opening hatch. A small 0.1 m diameter pressure relief vent was placed at the floor 


level to relieve the pressure from the initial burner ignition. Without this vent, the initial 
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pressure rise sufficed to activate the pressure relief panel. A computer controlled cover 
closed over this vent 15 s after ignition. Additional apparatus description can be found in 
refs. 5 and 7. 

6.2.2 Species Concentration: In order to characterize the compartment 
conditions prior to a backdraft, the species concentration histories in the upper and lower 
layers were recorded. Gas concentrations measured were: oxygen (O,), carbon dioxide 
(CO,), carbon monoxide (CO), and total hydrocarbons (HC). Continuous gas samples 
were taken with stainless steel probes. Only one analyzer of each type was available so 
the lower layer concentrations were measured in only 3 of the 11 experiments. For the 
upper layer concentration measurements the probes were located 0.6 m from the opening 
wall, 0.2 m from the ceiling, and 0.6 m from the side wall, as shown in Fig. 6.1. For the 
lower layer species measurements, the probes were located 0.15 m above the floor and 
0.76 m from the burner wall. The O,, CO, and CO, samples were taken through an 
unheated sample line in which the soot and water vapor were removed by glass fiber and 
desiccant filters, respectively. 

The HC concentration was more difficult to measure. The hydrocarbon sample 
required a separate heated sample line to prevent the loss of hydrocarbons due to 
condensation. A flame ionization detector was used to measure the hydrocarbon 
(methane) concentration. The effective range for this meter was 0 to 1%. The expected 
range of hydrocarbons was of the order of 25%, by volume. It was therefore, necessary to 
dilute the sample. The dilution system was designed to mix the compartment sample with 
heated ambient air in a ratio of 25 to 1. The flow rate of the dilution air and the sample 
were determined by measuring the pressure drop over a fixed length of tubing, 1.7 mm ID. 
Sample flow rates were monitored continuously during the experiment. Typical dilution 
ratios would vary from 22 to 25 depending on the compartment gas temperature and the 


soot build up in the sample line. A hot water jacket around the hydrocarbon sample line 
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and the dilution air line kept the gases over 60°C, well above the maximum calculated dew 
point of 44°C. To obtain the final species concentrations in the upper layer it was 
necessary to calculate the upper layer H,O concentration, since there was no direct 
measurement of H,O. Detailed species balances were performed’. 

6.2.3 Temperatures: A vertical thermocouple tree was placed 0.6 m from the 
opening wall and 0.2 m from the pressure relief panel, as shown in Fig. 6.1. The 
thermocouples were made from 0.5 mm type K thermocouple wire with a stainless steel 
overbraid. The average bead diameter was 1.1 mm. The ten thermocouples were located 
at 0.10 m intervals, with the highest thermocouple at 0.15 m below the ceiling. The 
temperatures reported here are uncorrected values. 

The thermal interface height history was calculated from the time dependent 
temperature profiles recorded from the thermocouple tree. The profiles were converted 
into unsteady idealized upper and lower layer temperatures using the method Quintiere et. 


al.10 applied to steady state temperature profiles: 


h, 


i(2 }ic= tm, Yr +, /T, 61) 


h 
and = {Tdx=[h, —h, JT +h, T, (6.2) 
0 


where T’ and T’” are the upper and lower layer temperatures, and h, and h,are the 
heights of the compartment and the layer interface, respectively. Equation (1) is a mass 
balance and Eq. (2) retains the same mean temperature as in the data. Assuming that the 
lower layer temperature was the arithmetic average of the two lowest thermocouples 
allowed the upper layer temperature and thermal interface location to be calculated from 


Eqs. (1 and 2). 
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6.2.4 Compartment Pressure: The compartment pressure history was recorded 
using an electronic pressure transducer with an effective range of 0 to 1250 Pa. The 
pressure port was mounted in the stationary wall opposite the pressure relief panel at floor 
level. The ambient pressure reference was taken outside the building. 

6.2.5 Opening Flow: The flow in and out of the compartment after the hatch was 
opened was recorded using six bidirectional probes in the opening. The probes were 
located in the horizontal center of the opening and 65 mm apart. The top and bottom 
probes were 43 mm from the soffit and sill. Probe velocities were calculated using the 


relationship given by McCaffrey and Heskestad!!: 


ORE) Pee | (63) 


where v is the velocity at the probe, Ap is the pressure difference measured across the 
probe, p is the density at the probe, and C(Re) is an empirical calibration constant which is 
a function of the Reynolds number. 

The temperature was recorded at each bidirectional probe in the opening using a 
bare bead 0.05 mm Type K thermocouple with high temperature glass insulation. Four 
aspirated thermocouples, designed according to Newman et. al.,!2 were installed to correct 
the temperatures for econ Data from the aspirated thermocouples revealed that no 
radiation correction was required. 

The total mass flow was calculated by integrating the velocity and density profiles 


over the height of the opening as: 


rh, = Ch*[p(y,t): v(y, t): bly, (6.4) 


tio =Ch*Lo(y,t)-v(y,t)-b]dy (6.5) 
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where m,,, and m,,, are the mass flow in the lower portion and the mass flow out the 
upper portion of the vent opening, C is the flow coefficient taken as 0.68, b is width of 
opening, and h,, h,, and h, are the height of the sill, soffit, and neutral axis, respectively. 

h, was determined empirically from the probe velocity data. The density was corrected for 
temperature using the thermocouple measurements. When flames reach the opening the 
thermocouples cannot respond fast enough to accurately measure the temperature and the 
mass flow calculations are terminated. 

6.2.6 Data Acquisition System: Data from each sensor were recorded using a 
HP VECTRA 80486-33 computer with a 8 channel multifunction analog and digital 
input/output board. Two 32 channel analog input muliplexors were connected to this 
system. A total of 31 thermocouple and 17 voltage channels were used. To increase the 
scan rate, data were written to a RAM drive and then downloaded to a file on the hard 
disk immediately following the experiment. The system was capable of recording each 
channel 50 times a second. For experiments greater than 600 s the data was collected at a 
rate of 10 scan/s until 20s before opening when the rate was automatically increased to 50 
scan/s. The reduced scan rate for the initial period of the experiment simply to reduce the 
data file size. Files for these experiments ranged from 4 to 10 Mbytes. 

In addition to recording the data, the computer also controlled the experimental 
procedures using solid state relays activated by a digital input/output board also installed 
on the computer bus. The computer controlled systems included: experimental clock, 
burner pilot light, fuel flow, vent cover, spark ignitors, opening hatch, and still camera. 
Each system had a manual override. | 

6.2.7 Procedures: Before each experiment a 60 s baseline was taken to record the 
initial conditions. A pilot flame was ignited at the burner 5 s before the start of the 
experiment. At 0 s a solenoid was opened on the methane flow to the burner and the 


clock was reset to zero. The burner was left on for a predetermined time period. Gas 
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flow to the burner was terminated 5 s before the hatch was opened. After opening there 
was a predetermined time delay ranging from 15-30 s before the spark ignitors were 


activated. 


6.3 Experimental Results 

6.3.1 Exemplar Data: The experimental parameters for the results presented in 
Figs. 6.2 - 6.5 are: gas flow rate of 70 kW, burner flow time of 535 s, hatch opening at 
540 s and backdraft ignition delay time of 20 s. Figure 6.2a shows the upper layer species 
mass fraction histories for O,, CO, CO,, and HC. Lower layer species concentrations are 
given in Fig. 6.2b. Idealized layer temperatures and height calculated from Eqs. (6.1 and 
6.2) are shown in Fig. 6.3a. For the first 70 s, the temperatures in the compartment rise as 
the layer within the compartment descends over the fire. The upper layer O, 
concentration is dropping as the fire consumes the available O, and the CO, concentration 
increases as a result of the combustion. Species concentrations in the lower layer are 
relatively unaffected during this early period. Over the time period of 70 s to 150 s the 
temperature in the compartment drops as the fire oscillates and then dies. After 150 s the 
fire is completely out and the compartment begins to cool as seen in the exponential decay 
in the temperature profile. Over this same period the upper layer HC concentration rises 
as the O, and CO, level off. The O, concentration in the lower layer drops considerably 
as the CO, and HC concentrations increase. After 240 s the hydrocarbon concentration 
in both the upper and lower layers is steadily increasing with the lower layer at ~20% that 
of the upper layer. The oxygen concentration slightly increases as air leaks into the 
compartment and the CO, concentration declines as compartment gases are lost due to 
leakage. At 480 s the HC analyzer became saturated and the slope of HC concentration 
approaches zero. A linear extrapolation would give a more accurate HC concentration in 


the last 60 s. At 540 s the hatch is opened and a gravity current enters the compartment. 
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The gravity current travels across the floor and is reflected back to the opening. After the 
gravity current reaches the opening oxygen rich air continues to pour into a new lower 
layer as the upper layer gases continue to exit the compartment. At 560 s the sparks are 
turned on and ignite a flame which travels through the compartment, driving the 
flammable gases out the opening and culminates in a large external fireball, ~5 m diameter. 
The spikes in the temperature shown in Fig. 6.3a are caused by the wave propagation 
through the compartment. Gas concentrations after 560 s are unreliable due to the highly 
transient effects of the backdraft. 

An idealized two zone approximation is compared with the temperature data from 
the thermocouple tree at 60 s and 540 s in Fig. 6.3b. At 60 s the compartment 
temperature is near the peak temperature and there is a substantial temperature gradient in 
the upper layer. The layer interface is located just above the burner surface. At opening, 
540s, the layer interface is still near the top of the burner and the compartment 
temperature is nearly uniform vertically indicating that a one zone approximation is 


reasonable at that time. 


106 


Mass Fraction 


Time (s) 


Figure 6.2a - Upper layer species concentration histories from backdraft experiment for O, 
( ), CO (---), CO, (------), and HC (~—--). Run 9 in Table 6.1 and 6.2. 


Mass Fraction 


Figure 6.2b - Lower layer species concentration histories from backdraft experiment for 
O, ¢ ), CO (-—-), CO, (------), and HC (——-). Run 9 in Table 6.1 and 6.2. 
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Figure 6.3a - Idealized two zone upper layer ( 
(.-----). Run 9 in Table 6.1 and 6.2. 
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Figure 6.3b - Idealized two zone approximation compared with temperature data at 60 s (- 
----,O) and 540 s ( ,0). Run 9 in Table 6.1 and 6.2. 
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Figure 6.4 shows the effects of a backdraft, i.e., short histories for the 
compartment pressure and the total mass flow into and out of the opening. The graph 
begins at the opening time, 540 s, and shows a 30 s period during which the gravity 
current enters the compartment, a backdraft is ignited and a large fire ball exits the 
compartment and conditions relax to quasisteady equilibrium. At 540 s the compartment 
pressure drops as the hatch falls open. There is a 20 s delay as the oxygen rich air enters 
the compartment. At 560 s the backdraft is ignited, arrow 1 in Fig. 6.4. Flames exit the 
compartment at 560.9 s as the pressure reaches 160 Pa, marked with arrow 2 in Fig. 6.4. 
The pressure spike in the compartment is caused by the large fireball erupting outside the 
compartment in the large, 900 m3, indoor facility with approximately 11 m? of vent area. 
The total mass flow into and out of the compartment is calculated from Eqs. (6.4 and 6.5). 
Negative values indicate flow into the compartment. After opening and prior to ignition, 
the total mass which flows into and out of the compartment steadily increases. After 
ignition, the flow into the compartment is reversed and all of the flow is out of the 
compartment as indicated by the period of zero slope on the inflow curve and the steep 
increase in the slope of the outflow curve, as shown in Fig. 6.4. Once the flame leaves 
the compartment and the fireball has subsided, the mass flow rate into the compartment 
reaches a quasisteady state as indicated by the constant slope in the last 8 s of Fig. 6.4, 


which is flow into the compartment as it cools. 
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Figure 6.4 - Compartment pressure ( ) and total mass inflow (------) and total mass 
outflow (——-) histories. Arrow 1 indicates ignition of the backdraft and arrow 2 


indicates flame out the opening. Run 9 in Table 6.1 and 6.2. 


Figure 6.5 is a series of video images from this experiment taken at 0.2 s intervals 
after ignition. The flame is initially spherical in shape and then stretches toward the 


opening, exiting the compartment 0.9 s after ignition. 
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6.3.2 Summary: Table 6.1 is a summary of the 11 experiments reported here. 
Column 1 is the experiment number. Columns 2 and 3 are the bummer characteristics, i.e., 
the burner flow rate and the time the burner gas is flowing. Columns 4 - 6 are results 
calculated from Eqs. (1 and 2) for the lower layer temperature, the upper layer 
temperature, and the layer height at opening. Column 7 is the observed ignition delay time 
between compartment opening and ignition of the backdraft. Column 8 is the peak 
pressure measured inside the compartment. Column 9 - 10 are the total mass which flows 
into and out of the compartment after hatch opening and prior to ignition, respectively. 
Column 11 - 12 are the total mass which flows into and out of the compartment after 
opening and prior to flames exiting the compartment. Column 13 is the run name which 


corresponds to the data given in appendix C. 
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Table 6.1 Summary of the 11 window experiments reporting the burner characteristics, layer temperatures 


and height, ignition delay, peak compartment pressure, and mass flow into and out of the compartment. 
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The burner flow rate was difficult to control and some fluctuation is seen in the 
data reported in column 2. Burn times ranged from 420 s to 540 s. Times greater than 
540 s were considered too hazardous to attempt safely. The ignition delay ranged from 15 
s to 30s. The position of the gravity current within the compartment can be calculated as | 
a function of time from salt water modeling assuming a constant Froude number discussed 
in chapter 3. Velocity current velocity for these experiments was approximately 0.5 m/s. 
For these experiments the gravity current reaches the rear wall approximately 7 s after 
Opening and returns to the front wall approximately 14 s after opening. The species 
concentrations are nearly constant for O,, CO, and CO, indicating the repeatability of the 
system. The HC concentrations histories are similar to Fig. 6.2a with a long gradual build 
up, although the slope of the curve changed slightly according to the burner flow rate. 

The idealized layer temperatures decreased as the burner times increased due to the energy 
loss to the boundaries and leakage of the hot layer gases. Layer interface height may be 
taken as approximately constant throughout the experiments. The maximum pressure and 
size of the fireball can be considered as a gross measurement of the intensity of the 


backdraft. They are seen to increase with the HC concentration as expected. 


6.3.3 Chemical Energy Accounting: An energy budget on the compartment was 
prepared to determine the amount of energy available for the backdraft. The integral over 
time of the burner flow rate can be calculated to determine the total chemical energy 


released into the system: 


Ex = |, (AHc- tig) dt (6.6) 


where E, is the total chemical energy released by the burner, to anytime t, AH. is the heat 


of combustion and mg, is the burner flow rate. The burner energy is assumed to be either 


consumed in the burner flame, stored within the compartment, or lost due to leakage in the 
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compartment boundary. The energy stored in the compartment at each time t is calculated 


assuming the compartment can be divided into two idealized layers as discussed above. 
Eta AL AON colle Goede), (6.7) 


where E, is the amount of chemical energy stored within the compartment, Yyo and Yyc 
are the mass fractions for the hydrocarbons in the upper layer and lower layer, 
respectively, and mY and m/ are the total mass of the upper and lower layers, 
respectively. The total mass in each layer is calculated from the layer height and average 


layer density: 
m™ =p” -w-L-(h,—h,), (6.8) 
m~=p"-w-L-h,, fort<to, (6.9) 


where ty, is the time when the compartment hatch is opened, p and p"" are the densities 
of the upper and lower layers, respectively, w is the compartment width, L is the 
compartment length, and h, and h, are the compartment height and the height of the layer 
above the floor. The energy lost by leakage is found by the difference, E, =E,-E,-E,, 
where E, is the energy lost from leakage and E, is the energy combusted at the burner. 


The leakage mass loss rate m, is then estimated from: 


E, = | AH. Yyernydt. (6.10) 


The chemical energy histories for the compartment are shown in Fig. 6.6 for run 9 


in Table 6.1 and 6.2. The dashed line indicates the total chemical energy which enters the 
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compartment through the burner, (Eq. 6.6). The solid line indicates the chemical energy 
which is stored within the compartment as a function of time. The chemical energy 
supplied by the burner simply increases linearly due to the constant flow rate (74 kW) 
through the burner. At 535 s the burner is shut off and the total chemical energy supplied 
by the burner remains constant. The total stored chemical energy is initially zero since all 
of the chemical energy released by the burner is consumed by the fire. As the fire burns 
out, the hydrocarbons start to build as unburned fuel is stored within the compartment. 
The stored energy continues to increase to 480 s. The difference between the burner 
supplied energy and the stored energy is a function of the chemical energy consumed by 
the fire and the chemical energy lost due to leakage out of the compartment. The amount 
of chemical energy consumed by the fire is approximately Ec~7 MJ. The chemical energy 
lost due to leakage is estimated at 9 MJ. From Eq. (6.10) that suggests a constant leakage 
rate of 5 gm/s. The total chemical energy history with that mass loss rate is shown by the 
dotted line in Fig. 6.6. With E,~39 MJ the energy available at the time of opening is E,~ 
19 MJ. 
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Figure 6.6 - Chemical energy histories for the energy within the compartment (. -) and 
the chemical energy supplied by the burner (— — —) and the total energy history with an 
assumed constant mass leakage rate of 5 gm/s (------). From run 9 in Table 6.1 and 6.2. 
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The available energy for backdraft is assumed to consist of: 1) the amount of 
energy which leaves the compartment prior to ignition, E,, 2) the amount which leaves the 
compartment after ignition and before flame exits the opening, E,, 3) the amount which is 
contained in the compartment when the flames reach the opening and not consumed there, 
Ey. Each of these will now be estimated. 

After opening the compartment but prior to ignition of the backdraft, the mass 
flow out of the compartment is assumed to be out of the upper layer and the flow into the 
compartment is assumed to be into the lower layer. The unburned chemical energy that 


has flowed out of the compartment is: 


E,=AHo. Yuo C bf" ['ply,t) v(y,t) dy dt, forto— t= t.(6.11) 


where Yio is assumed constant over that time interval. 
After ignition of the backdraft, the compartment gases are assumed to be well 
mixed and the hydrocarbon mass fraction is taken as: 


Dl UL Lb ALL 
Yuc (Met, — Mout) + Yue * Mat, 


UL LL 


Vio 
My=t, + Mat, + Min — Mout 


for t;> t > to, (6.12) 


where t,, is the time when flames exit the compartment, Y;,¥0*, is the mixture hydrocarbon 
mass fraction. Therefore, after ignition and before flames exit the compartment, the 


energy exiting the compartment ts: 


E,=YM* Cbf'* ['[o(y,t)- v(y,t)] dy dt, for t;>t> ta, (6.13) 


To estimate the amount of chemical energy within the compartment which is not 


consumed it is assumed that all the oxygen within the compartment when flames reach the 


117 


opening is burned. The amount of oxygen which is stored within the compartment at 


opening, mg? , 1S, 
03, 4/5,UL ee Eben 
ms? =(Yo, -m'~+Yo,-m~), (6.14) 


where Mon and Gr are the mass fraction of oxygen in the upper and lower layer 


respectively. After opening, oxygen is gained from the incoming air and mass 1s lost in the 


out flow from the upper layer for a net change in the mass of oxygen, Am?? , given by 


. phy 
Amy? = Yeas bf if p(y,t) v(y,t) dy dt— Ygne bf if p(y,t) v(y,t)dydt (6.15) 


for tp > t > t,, and where See is the oxygen mass fraction of air. 


After ignition, the compartment gases are assumed to be well mixed as in Eq. 
(6.12) and the oxygen mass fraction is taken as: 


UL UL LU LL 
Yo, (Meat, — Mou) + Yo, Midi, 


UL LL 2 


You = 
Mi=t, te My=t, + Min — Mout 


nee ee th, (6. 16) 


where YOu is assumed constant over that time interval. The change in the stored oxygen 


concentration over the interval t;> t > tg, Am>3, is then taken as: 


i tro (Ba ix ro (B, 
Amps = Y3C bf” [)"ply,t) vy, t) dy dt— YORC bf” p(y.) v(y,t) dy dt (6.17) 


The amount of oxygen that is stored in the compartment when flames exit the 


compartment is then, my3, then: 


Oma to ) O 
m6 = my? +Am,? +Am;6, (6.18) 
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The maximum amount of stored chemical energy at the time flames reach the 


opening which can burn in the compartment, E,, is, 


Bo = les meas (6.19) 


where the 12.5 factor is the amount of energy released per unit mass of oxygen based on 


an oxygen-methane reaction. The energy is the compartment, but not consumed in, 


The portion of the chemical energy stored in the compartment at opening that 


remains available for the fireball, E,,, can now be estimated as: 


Epp = E; +E; +Ey = Es —-Eo, (6.21) 


In Fig. 6.6 for run 9 in Tables 6.1 and 6.2 at 540 s the hatch covering the opening 
falls open and the upper layer gases start to exit the compartment. As these gases leave 
the compartment the stored energy, E, = 19.4 MJ, drops as demonstrated by the steep 
negative slope seen in Fig. 6.6 after 540 s. At ignition, t= 560 s, E, = 5.1 MJ and the 
slope of the stored chemical energy line becomes very large as the compartment gases are 
forced out of the compartment by the expanding flame, E, ~ 2.7 MJ. Some of the 
remaining chemical energy stored within the compartment has been consumed there E, ~ 
5.3 MJ; Ex, ~ 6.3 MJ was not. The chemical energy forced out of the compartment, E,, ~ 
14.1 MJ, may burn as a large fireball outside of the compartment. The combustion 
efficiency of the fireball remains unknown. Table 6.2 is a summary of the calculated 


chemical energy parameters. Column | is the experiment run number. Columns 2 - 5 are 
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the compartment species concentrations at opening for O,, CO, CO,, and HC, 
respectively. Species concentrations printed in parenthesis indicate the sample was taken 
from the lower layer. The average hydrocarbon mass fraction for the compartment 
evaluated at backdraft ignition is shown in column 6. Column 7 is the chemical energies 
stored in the compartment when the compartment was opened. Column 8 is the chemical 
energy which flows out of the compartment prior to backdraft ignition. The chemical 
energy which flows out of the opening between compartment ignition and the flames 
exiting the opening is shown in column 9. Column 10 is the amount of energy bumed in 
the compartment assuming all the available oxygen is consumed. Column 11 is the energy 
in the compartment which is not burned there. Column 12 is the amount of chemical 
energy which is available for the fireball. Column 13 1s an estimate of the diameter of the 
fireball which burns outside the compartment. Column 14 is the run name which 
corresponds to the data given in appendix C. 

The presence of a fireball is used as the indication of a backdraft in these 
experiments. When the fireball is not observed, the flame travel in the compartment is 
noticeably slower and less intense than in the experiments with a fireball. The fireball also 
gives a crude indication of the severity of the backdraft. If the fireball were contained 
within a volume scaled with the compartment, such as a corridor, the pressure rise within 
the compartment would be significantly higher, possibly enough to cause structural 


damage. 
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6.4 Future Research 

Future work should concentrate on designing and building a full scale apparatus 
using a controlled fire source and possibly more realistic fuels. Improved gas analysis 
would also be useful to determine what hydrocarbons are present in the real fuel 
experiments. Additional opening geometries using different vent/wall area ratios and 
locations could also be used to investigate the effects vent flow has on the backdraft 


severity and chance of occurrence. 
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CHAPTER 7 
CONCLUSIONS 


7.1 Summary 

The backdraft scenario presented here accurately describes the phenomena known 
as backdraft. Experiments show that the ignition of the backdraft does not occur 
immediately upon opening the compartment. It is this time delay caused by the gravity 
current propagation which creates a hazard to firefighters who enter a compartment and 
become trapped in the backdraft process. 

The gravity current entering a compartment is both qualitatively and quantitatively 
similar to other naturally occurring gravity currents. The nondimensional velocity for the 
full opening compares well with the perfect fluid theory. The structure of the entering 
gravity current head for the full open condition shows a shallow mixed region nding in the 
current near the interface. It also has the detailed features reported for steady state gravity 
currents, i.e. billows, lobes, and clefts. Entering currents for the slot, window, and door 
openings show a different gravity current structure with the mixed region occupying 
nearly the entire current due to the enhanced mixing near the opening. Similar detailed 
features of the steady state gravity current appear in these currents. 

The values of v* and h* obtained here for a variety of opening geometries, are 
independent of the density difference ratio, 8. The exiting gravity current is also 
independent of B, and has a v* approximately equal to the entering current. These results 
can be applied to predict the time to ignition for a backdraft with compartment and 


opening geometries similar to the conditions reported here. 
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The two dimensional numerical simulation performed at the National Institute of 
Standards and Technology accurately predicts the compartment gravity current for the 
fully open condition. The two dimensional computational density profile shows that the 
main body of the current is made up ambient fluid while a mixed layer made up of large 
vortical structures travels along the shear interface between the two fluids. Such structure 
is also observed in the salt water experiments although the individual vortices are not 
visible due to the averaging inherent in of the flow visualization techniques. For the slot 
opening, the results of the two dimensional simulation also compare favorably with the 
experimental results. The computed density profile shows that the structure of the gravity 
current for the h,/3 centered slot opening 1s significantly different from the fully open 
condition. For the slot opening, the entire gravity current is filled with a complex 
structure of large scale vortices. This large scale mixing for the h,/3 centered slot opening 
results from the rearward facing step, forced by the lower edge of the slot. Two 
dimensional transit times predicted by the numerical simulation compare well with the salt 
water experiments although the computed times are slightly longer for the slot opening. 
Comparing the computational transit time with the time to ignition for the backdraft 
experiments gravity current “iba shows good agreement. Numerical velocity profile 
predictions in the opening compared well with backdraft experiment results and fall well 
within the experimental error bounds. 

The results presented for the quantitative backdraft experiments show that the 
hydrocarbon concentration must be >10% in order for a backdraft to occur. When the 
hydrocarbon concentration is < 10% the flame travel is slow and the compartment 
overpressure is much lower. As the hydrocarbon concentration increases, the 
compartment overpressure increases and the backdraft becomes more severe. 


Hydrocarbon concentrations > 15% exhibit large fire balls outside the compartment and 
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would result in significantly higher pressures if contained within an adjacent corridor or 


room with a volume of the order of the experimental compartment. 
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7.2 Future Research 

The salt water modeling was limited to only one compartment aspect ratio and 
limited opening geometries. Future work should focus on different compartment aspect 
ratios, smaller openings, and openings at the floor level. Improved flow visualization 
techniques including slit lighting to obtain a two dimensional image of the gravity current 
are recommended. Time dependent concentrations should be measured within the gravity 
current for comparison with computational results. 

Future backdraft experiments should concentrate on two areas. The first should 
focus on designing and building a full scale apparatus i.e. a minimum h, = 2.4 m. Full scale 
experiments should focus on quantifying the effects of backdraft and how to prevent or 
reduced the impact of backdraft. Suggested areas to investigate include: 1) Can vertical 
ventilation prevent backdrafts? 2) Effects of horizontal ventilation when it is not possible 
to vent the ceiling of the structure. 3) Impact of closing the door after entering a 
compartment. 4) Effect of containing the fireball in an adjacent space of the same order as 
the backdraft compartment. 

Further research is also necessary on common fuels which may contribute to 
backdraft. Results presented here indicate that high unburned fuel concentrations are 
necessary in order for a backdraft to occur. Experimental work is required to determine 
the effects of ventilation on the pyrolysis rate of the fuel and to quantify the amount of 
pyrolyzates consumed in the flame as a function of the available oxygen concentration. 
Gas analysis capable of determining the nature of the hydrocarbons present would also be 


helpful in calculating flammability limits. 
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APPENDIX A 


SPECIES CONCENTRATIONS 


In order to evaluate the composition of the upper layer of the compartment it is 
necessary to determine all of the major species present, specifically O,, CO,, CO, HC, 
H,O, and N,. Direct measurement of species concentrations were only possible for only 
O,, CO,, CO, and HC. Concentrations of H,O and N, had to be calculated using species 
balances on O, and total C. The following assumptions were necessary for the analysis: 1) 
that the upper layer is well stirred and 2) that the overall reaction was: 

CH, +a(O, +3.77N,)+bH,0 > cCO, +dCO0+eH,0+fN,+gCH, (A-1) 
Using these assumptions, which are reasonable since the fire source was a gas burner 
flowing methane, an overall balance on the oxygen and total carbon will yield the two 
concentrations, H,O and N,. 

The first step is to determine the amount of N, which is present in the system. 
This is done by a simple O, balance. Starting with the dry sample which has had all of the 


water scrubbed out using a desiccant filter the total amount of oxygen is calculated: 


ne =n“ (Xe + ies +2X') (A-2) 
Symbols follow standard chemical notation when possible and a complete listing appears 
in the nomenclature. The oxygen is not only in the analyzed gases but also in the water 
which has been scrubbed. The amount of H,O which is present is a result of the 


combustion of CH, and can be estimated from the following relationship derived from a 


C balance, 


Hae = ZT =f n“co) — Te Xen =f ee (A-3) 


128 


Combining Eqs. (A-2 and A-3), multiplying by 3.77, i.e. the molar ratio of N,/O., and 


assuming np=1 gives the total amount of N,: 


gs, Te See + X55] (A-4) 


Knowing the amount of N,, the remaining unknown is the concentration of H,O. 
There are two sources of H,O, that produced by combustion and the amount present in 
the air. The amount of H,0 produced by combustion has already been estimated in eq. 
(A-3). The amount of H,O present in the air is obtained from the following definition of 


relative humidity: 


we fo 622 ai eo (A-5) 
(101, 325- OPsart) 


where w is the specific humidity, is the relative humidity, and p,,, is the saturation 
pressure (Pa). The relative humidity was recorded before each experiment and p,,, can be 


calculated from the following relationship taken from ref. 1: 


ft —-A 
Pedi) = Cee eZ | > (A-6) 
Vv 


where A=3.18 x 10° J kg!) B=2470 J kg! KK", and C=6.05 x 1026 N m?. Knowing the 
relative humidity and temperature the amount of H,O in the ambient air can be calculated 


from the following relationship: 


Seen (A-7) 


n any 
H,0 
Mite 


The mass of the dry air is calculated from the: 
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Mee n*| Mg [xs + 2X00: 3x, 4M,.x5, | (A-8) 


Using eqs. (A-3 and A-7) and the assumption that np>=1, the total number of moles in the 


system can be calculated as: 


eae, c Air Cc Air 
no=n +n +n, ,o— itn tn, (A-9) 


Applying eq. (A-9) the mole fractions for the upper layer can be calculated as, 


AxrA 
x x n > Fp 
oO, EL 
n 
AxzA 
xe e 10 
co, nt 
AyzA 
x n Xo 
CORm i 
n 
A 
aed Se 
Now aa n° 
xt ieee 
HO nt 


The final result are reported as mass fractions. To convert the mole fractions the 


molecular weight of the upper layer sample is calculated, 


M*= XN a Kee ie foNes Mears 2S dN OG Ti Spat : (A-10) 


Using the definition of mass fraction we can convert mole fractions to mass fractions with: 


ait (A-11) 
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APPENDIX B 


COMPARTMENT DETAILS 
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Pressure Relief Panel Detail 


1.2 mm (18 gauge) 
0.1 m steel studs 
nominal spacing 0.46 m 


galvanized steel sheet 
0.047 mm thick (18 gauge) 


one 25 mm (1”) thick layers 
ceramic fiber insulation 
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Wall Detail 


1.2 mm (18 gauge) 
0.1 m steel studs 
nominal spacing 0.4 m 


two layers 16mm thick (5/8") 
type X fire rated 
gypsum wallboard 


RCRA RS 


two 25 mm (1”) thick layers 
ceramic fiber insulation 


Ceiling Detail 


1.2 mm (18 gauge) 
4” steel studs 
nominal spacing 0.46 m 


two layers 16mm thick (5/8”") 
type X fire rated 
gypsum wallboard 


<i 


ELIE LEE LEE 
RRA RAGS 


two 25 mm (1”) thick layers 
ceramic fiber insulation 
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APPENDIX C 


BACKDRAFT EXPERIMENTAL DATA 
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Two Zone Approximation -Temperature & Interface Histories 
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Total Mass Out (kg) 


Layer Temperature History - TC Tree Data 
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Two Zone Approximation -Temperature & Interface Histories 


850 
800 
750 
700 
650 
600 
550 
500 
450 
400 
350 
300 = 
250 


Temperature (K) 


0 60 120 


Time (s) 


Layer Height (m) 


ile atts TOs oe cever 


Mass Fraction History - Compartment Gases 


Mass Fraction 


0 60 120 


Time (s) 


YO 2 eV (COlm ars se aie 
Run: P3EXP27 


140 


Y(HC) 


Bidirectional Probe Velocity History 


Velocity (m/s) 


Pressure (Pa) 
Total Mass Out (kg) 


Sevaleera 
a), 


ata am BGA noble tate Data 


180 185 190 195 200 205 210 
Time (s) 


Comps Pe". _. .eMass’ Out (o> = Ath Mass In 


Run: P3EXP27 


14] 
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Two Zone Approximation -Temperature & Interface Histories 
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Bidirectional Probe Velocity History 
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Two Zone Approximation -Temperature & Interface Histories 
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Bidirectional Probe Velocity History 
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Two Zone Approximation -Temperature & Interface Histories 
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Bidirectional Probe Velocity History 
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